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ON SOME PROPERTIES OF GENERAL MANIFOLDS RELATING 
TO EINSTEIN’S THEORY OF GRAVITATION. 


By J. A. SCHOUTEN AND D. J. StrrRuIk. 


In two very interesting communications in this JouRNAL* Mr. E. 
Kasner recently deduced some properties of the four-dimensional manifolds 
for which G,,, = 0. Now it is possible to demonstrate two of these proper- 
ties in a very short way for n-dimensional manifolds without calculating 
any three-index-symbol, as a direct consequence of some general theorems. 
Further it is very easy to give a generalization for manifolds, which obey 
not Einstein’s gravitational equations but the more general equations of 
de Sitter for a quasi-spherical world. 

By V,, is meant a general manifold of n dimensions, by R, a euclidean 
one, by S, a manifold with constant Riemannian curvature and by C, a 
manifold which is conformally representable on an R,. 

1. In a recent publicationt one of us showed that the necessary and 
sufficient condition for n > 3, that a V, isa Cn, is that the curvature-tensor 
(Riemann-Christoffel-Tensor) B,,,, can be written in the form: 


1 


P, = Ay, An 
where L,, is any symmetrical tensor. From this we conclude 


1 Loo”, 


Lp = — Cw + 
Hence B,,,« vanishes for a C,, if G,, vanishes: 
Of all C,, the only ones for which G,, = 0 are Rp. 
For n = 3 B,,,. has always the form (1), hence: 
Every V3 for which G,, = 0 is an R3. 
This theorem can be generalized in the following way. If 
(4) Gy = 


* AM. JOURNAL OF MATHEMATICS, 43 (1921). 

+ J. A. Schouten, “Ueber die konforme Abbildung n-dimensionaler Mannigfaltigkeiten 
mit kwadratischer Maszbestimmung auf eine Mannigfaltigkeit mit euklidischer Masz- 
bestimmung,” Math. Zeitschr., 1921, pp. 58-88. 
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214 SCHOUTEN AND StRUIK: Properties of General Manifolds. 


(in relativity the equation of de Sitter for a quasi-spherical world) from (3) 


follows 
©) 2(n — 
and thus: 


2(n— 1) 


Now it is a well-known theorem that B,,,, being of this form the Vy is 
an S,: 
Of all C,, the only ones for which G,, = Agu, are Sn. 
Every V3 for which G,, = dg,, is an 83. 
2. If a V, is imbedded in a Vy41, and if the indices yu, v, p, belong to 
the fundamental variables in the first manifold, then the generalized 
theorem of Gauss is: 


where B,,,. is the curvature-tensor of the Vn, the curvature-tensor 
of the V4; and h,, the (symmetrical) second fundamental tensor of the 
V,, with respect to the Vi4i.* Now if the is an Rayi, B’ = and 
from (7) we get: 


(8) Gu = 9 — h,,h, 
where 
(9) h = 


Now we introduce a system of rectangular coérdinates in the principal 
directions of the symmetrical tensor h,,, which are, according to (8), 
also the principal directions of G,,. Let the orthogonal components be 
denoted by the indices 2, 7 = 1, ---,n. Then G,, has n (real) components 
G;;, and h,, has n (real) components h;;, the components G;;, hij, t+ 7 
being zero, and (8) becomes 


(10) Gi; = hihi; hi; h;; 
j 
where no summations are to be taken that are not indicated by the sign 2. 
If now G,, = 0, we get from (10) that either all components h;; vanish or 
all components h;; but one, in both cases from (7) follows that B,,,. vanishes. 
Thus: 


*See the above cited paper p. 24, or Bianchi-Lukat, Vorlesungeniiber Differential- 
geometrie, first edition, p. 623, or G. Ricci, Atti Acad. Lincei (1902), p. 359. 
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3) No V,, with G,, = 0, that is not an R,, can be regarded as imbedded in an | 
Riu. The R, rs either flat or developable. | 
This theorem can also be generalized. If 


where G,,’ belongs to the Vn41, then (10) gets the form: 


A- = hihi; hi; 2 h;;. 
Jj 


For n = 2 hy and hy are not determined by these equations. For n + 2 


is it is easy to prove that the only possible solutions are: 
1 
hy = = hmm = ( )s 
| [m—-1, 
d m 
m=n— Mm; m = 0, 2, 4 for n = 4, 
m = 0, 2,3, °+:,n—2,nforn > 4. 
For these manifolds the generalized Codazzi-equationst give the follow- | 
‘ ing results: 
If for a Vin in hy = +++ = hmm ts a complete group of equal principal 
orthogonal components of the second fundamental tensor, then the corresponding 
congruences of principal curvature form a Vm with only umbilical points with 
respect to the If the Cnii is the value of = = hmm on 
each of the ©”-"V», remains constant. If hy, = +++ = mm remains constant i 
in the Vn, the o"-"V» are geodesic in the Vn. 
In consequence of these theorems and the above-mentioned solutions 


we can prove the theorem: 

A V, with G,, = dg, that is imbedded in an Snii, n> 3, ts erther a 
manifold with only umbilical points (a hypersphere) or it contains ©"—"V,, 
and also ©"™V ym with only umbilical points with respect to the Sn41, where ; 
m may have the values 2, ---,n— 2. The Vm are geodesic with respect to 
the V,, and the value of the principal radii of curvature remains constant on 
the Vn. 

Also it can be proved that each of the above-mentioned V’,, is imbedded 
in an S,,41 that is geodesic in the S,,,, and each of the V,» in a geodesic 


Hence a Vin S; with G,, = dg,, contains two systems of S2 with only 
umbilical points, which have all the same curvature. 


* For a V2 in S; these equations are identities. 
Compare Ricci, l.c. 
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3. The first theorem mentioned under 1 is a special case of a more general 
theorem.* 

If a V, with fundamental tensor g,, is represented conformally on a V, 
with fundamental tensor ’9,,: 


(11) "Vu i 
then we can deduce for the curvature-tensors B,,,, and 'B,,,. of both Vy: 
where 
(13) Pur = + — 28,819’ Juo 

From (11) follows: 
(14) "Guy = + — — NDyy)- 
Now when G,, = 0 and ’G,, = 0 from (14) we get p,, = 0, hence 
(15) 


When two V, with G,, = 0 are conformally representable on each other the 
Riemann-Christoffel curvature-tensors are equal. 
Detrt, Houuanp. 


*Schouten, ab. cit. p. 79. 
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GEOMETRICAL THEOREMS ON EINSTEIN’S COSMOLOGICAL 
EQUATIONS. 


By Epwarp KASNER. 


I wish to generalize here some of my results published in the January. 
and April numbers of the AMERICAN JOURNAL OF MATHEMATICs (Vol. 48, 


1921, pp. 20, 126), relating to Einstein’s original equations of gravitation 
(in space free from matter), 
(1) G., = 0. 


Later Einstein introduced a so-called cosmological term involving a con- 
stant , the equations being then 


(2) Gur — = O 
More recently* he has employed the form 
(3) = 0, 


where G is the scalar curvature g*°G,,. I shall refer to (3) simply as the 
cosmological equations. Every solution of the former equations (1) is, of 
course, a solution of the latter (3), but not vice-versa. The ten equations 
(3), as Einstein shows, involve one extra dependence as compared with the 
ten equations (1). 

§ 1. Five Dimensions. 


I shall first take up the question of dimensionality (that is, class of the 
quadratic form). In the April paperf it was shown that no solution of 
(1) can represent a 4-spread imbedded in a 5-flat (except in the trivial case 
where the 4-spread is euclidean, that is, has zero Riemann curvature). We 
now inquire what solutions of (3) can be imbedded in a 5-flat, and shall 
find that there are actually two distinct possibilities. 

Using the notation of the April paper, we write our spread in the form 


w = f(x1, Xe, 23, 24). 


Referring to the formulas on p. 128, we have, in the standard coérdinates 
there employed, involving the four principal curvatures k;, 


Giz = 0, ete.; Gu = ky (ke +hk3+ ks), etc. 


* Berichte Berlin Akad. d. Wiss. (1919). The complete equations when matter is 
present of course involve the energy tensor 7,,,. See also Kopff, Grundziige (1921), pp. 
163-165, where the author refers to the field equations of the first, second, third kinds. 

t “The Impossibility of Einstein Fields Immersed in Flat Space of Five Dimensions,’’ 
Vol. 48, pp. 126-129. 
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Also by the footnote on the same page, we have 
G=-— 2 (hike + + + + + 


Substituting in the cosmological equations and simplifying, we have the 
following set of equations for the determination of the four principal curva- 
tures, 


(4) ha) = + + = + + hee) = + ho + hs). 
Subtracting say the second member from the first, we have 
(4’) (ky — ke) (kz + ks) = 0, 


with five similar equations obtained by permuting the subscripts. Hence 


either 
(4’’) hy = hy or ks == hea, ete. 


If the common value of the four expressions in (4) is zero, then we have 
exactly the system (12) of the April paper, giving merely the trivial case 
where 3 or 4 of the k’s vanish (which means that the manifold is euclidean). 
Otherwise we find from (4’’) two possible types of solution 


(a) ky $0, 
(b) = ko = ks ke, (0). 


In the first case (a) the four principal curvatures are equal at every 
point, that is, every point is umbilical. It follows then from known 
theorems that the 4-spread must be a hypersphere. This, of course, 
checks up since it is known that a 4-dimensional hypersphere is actually 
imbedded in a 5-flat and is actually a solution of the cosmological equations. 
(It is sometimes referred to as DeSitter’s “Spherical World.’’) 

In the second possibility (6), we have, at every point, the four principal 
curvatures numerically equal, but two of them are positive and two are 
negative. We may say then that every point is semi-wmbilical. It is 
immediately seen that the Riemann curvature of such a spread is not con- 
stant; for we find, in our special coédrdinates, that the conditions for constant 
Riemann curvature are 


(5) kik kiks koks = — kskg. 


For the spherical solutions (a), these products are all equal; but in the new 
case (b), two of the products are positive and four are negative. We may 
term a spread of this new type, a hyperminimal spread, since we may think 
of it as a generalization of ordinary minimal surfaces, which have the 
property that the two principal curvatures are numerically equal, but 


opposite in sign. The actual existence of such hyperminimal spreads 
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depends on the consistency of a certain set of three partial differential 
equations, of the second order, which can readily be written down. It may 
be that no solution exists, or that all the solutions are imaginary, but the 
possibility still remains open. 

TuEoREM IJ. If a four-spread imbedded in a five-flat is to obey Exnstein’s 
cosmological equations (3), then either every point 1s umbilical (giving a hyper- 
sphere), or else every point is semi-umbilical (giving a possible type of hyper- 
minimal spread). 

§ 2. CoNFORMAL REPRESENTATION. 


I next take up the generalization of a theorem given in the January 
paper.* It was there shown that the only solutions of Einstein’s original 
equations (1), which have the same light equation as the euclidean or Min- 
kowski spread, are themselves euclidean; that is, if the spread is to be con- 
formally representable on a 4-flat, it must have zero Riemann curvature. 
We shall now prove: 

THEOREM II. The only spreads which obey the cosmological equations 
(3) and which are conformally representable on a 4-flat are those that have 
constant Riemann curvature (that is, the spread must be of spherical or pseudo- 
spherical character). 

For this purpose we use the notation of the earlier paper, in particular 
the formulas on pp. 22-23. Using the value of the tensor G,, there calcu- 
lated, we find that the scalar curvature is 


{INi;+ 2N?7}, where A= 
Substituting then in the equations (3) we find the following system of equa- 


tions for the determination of the unknown function NV: 
Nw = ete. ; 


3Ni1 — Noo — N33 — Nag — 3N2 + N2+ + NZ = 0, ete. 
Employing the same transformation N = — log M as in the earlier paper, 
we obtain this simple system: 
Mi. = 0, = 0, etc.; 
My, = Mo = M33 = Mas. 


The general solution is obviously 


(6) M = + x2? + 437 + + + + + + as. 


* “Hinstein’s Theory of Gravitation: Determination of the Field by Light Signals,”’ 
Vol. 48, pp. 20-28. Apparently without knowledge of this paper, Ogura has recently given 
this special theorem in Comptes: Rendus, Nov. 7, 1921. 

}In the January paper, by a.typographical error, an extra member }M2M;,? was 
omitted in the second line of the corresponding set (10), p. 23, but the final result there 
given is correct. 
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Our differential form is therefore, since \ = M~, 
(7) ds? = M~ (dx? + day? dx3? + dx,”). 


This is recognized as a manifold of constant Riemann curvature. (The 
curvature becomes zero only when the relation a;?+-a,’+a;°+a/?—4aa;=0 
holds, verifying the special theorem of the earlier paper.) 

Since a hypersphere can be mapped conformally on a 4-flat, it follows 
directly from Theorem II that the only cosmological solutions which can 
be represented conformally on a hypersphere are those of constant Riemann 
curvature. 

We may also generalize the discussion of approximately-euclidean mani- 
folds, given on pp. 27, 28 of the January paper, to approximately-spherical 
manifolds. The final result is 

THEOREM III. If two approximately-spherical spreads, both obeying the 
cosmological equations (3), admit conformal representation upon each other 
(thus having the same light equation), then they are necessarily isometric, except 
for a homothetie transformation. 


§ 3. SoLuTIoNS DEPENDING ON ONE VARIABLE. 


A simple example of a solution of (1) where the potentials involve only 
one of the variables is 
(8) ag — + dao? + dz’). 


All orthogonal solutions of this type (see Bull. Amer. Math. Soc., vol. 
27 (1920), p. 62) are easily shown to be reducible to the form 


ds? = + + + 
(9) dg +a3 +a,=a, +1, + a3? + ag = (a, + 1)?. 


This can be put in the static form, and is completely determined by its 
light rays. 

Analogous solutions of (3) are stated in an abstract printed in Science 
referred to below. We shall state the general result as 

THEOREM IV. All cosmological solutions for which the four potentials in 
the orthogonal form are functions of one of the four codrdinates can be found 
explicitly by elementary algebraic and transcendental functions. The cor- 
responding spreads can be imbedded in a 7-flat. 


§ 4. An ALGEBRAIC SOLUTION. 


We also state, omitting the easy proof, 
THEOREM V. [f the quaternary form ds* = g,,dx,dx, 1s to be expressible 
as the sum of two binary forms, one involving say x1, 2X2, the other involving 
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say x3, x4, and 2f the cosmological equations (3) are to be obeyed, then the only 
solution (except for a constant factor) is 


ds? = + dae”) + + 


This can be imbedded in a 6-flat with cartesian codrdinates (X,1X2X3X4X5X6), 
the finite representation being 


Ar Xo + — + Xs" ae =]. 


Excluding the obvious flat and spherical solutions, this is apparently 
the simplest solution of Einstein’s equations which has thus far been 
obtained, and is the first case where the finite solution is an algebraic 
spread. In the example (8) given in § 3, the potentials g,, are algebraic, 
but not necessarily the corresponding finite spread. 


CoLUMBIA UNIVERSITY, 
NEw YorK. 

* The theorems of the present paper were first published in Science, Vol. 54 (Sept. 30, 
1921), pp. 304-305. A typographical error on page 305 should be corrected as in formula 
(8) above. See also a paper appearing in the Mathematischen Annalen, entitled “The Solar 
Gravitational Field Completely Determined by its Light Rays.” An independent proof 
of Theorem II, together with very elegant proofs of two of my previous results, apply- 
ing to forms in 7 variables, is given in the paper by Prof. J. A. Schouten and Dr. D. J. 
Struick page 213 of this volume of this Journat. The authors were kind enough to 
send me a copy of their manuscript. 
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ON THE FERMAT AND HESSIAN POINTS FOR THE NON- 


The euclidean triangle has two Fermat or equiangular points. A Fermat 
point F is defined projectively by the condition that the lines joining F to 
the absolute point pair J, J are the hessian of the triad of lines joining F 


EUCLIDEAN TRIANGLE AND THEIR ANALOGUES 
FOR THE TETRAHEDRON. 


By C. M. Sparrow. 


to the vertices. If we take actual perpendiculars on the sides as coérdinates 
21, Xe, 23, the points F are transformed by the Desargues transformation 
x = 1/y into a pair of points H, which are defined most simply by the con- 
dition that the feet of the perpendiculars from H on the sides form an equi- 
lateral triangle. The Fermat points may be defined in another way which 
will be found very suggestive. If we consider a particle acted on by three 


equal forces directed along lines through the vertices, the points F will be 
positions of equilibrium. We may substitute geometrical for dynamical 
ideas by considering the potential energy. If 11, r2, 73 are the distances to 
the vertices the points F give the stationary values of 7; + r2- 73. The 
choices of sign, four in number, correspond to the different combinations 


of pushes and pulls, only two of the four having solutions. 
In the non-euclidean plane our point pair J, J is replaced by the conic 
(line codrdinates) 


where the c’s are the cosines of the internal angles, which in this case are 
Independent. We will denote by Ag the discriminant of ©. 
Lines from a point 24, 22, x3 to the vertices are 


and the hessian pair of this triad is (1/21, w/22, w?/a3) and (1/2x1, w?/a2, w/23) 
where w is a complex cube root of unity. If x is a Fermat point, this pair 
touches 2. The two equations thus obtained reduce to 


Writing y = 1/2 we get the pencil of conics 


The four base points of this pencil satisfy the definition of the hessian 


= + + £3? — — — = 0, 


1 1 2¢3 1 


2c 1 


(1) 
v3 


+ + 2c3y1Y2 yv + + 2coy1Y3 = + + 2C1Y2Y3- (2) 
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points H. We thus see that in the non-euclidean plane there are four 
Fermat points and four hessian points. Noting that the equation 
ao? + a3? + 2ci%273 = 0 gives the tangents from 1, 0, 0 to 2 and that each 
of the three conics is on a vertex, we see that each conic is on the vertices 
of a quadrilateral formed by two pairs of tangents and on the remaining 
vertex. If Ag = 0 two of the points H coincide with J and J and the conics 
are all circles. The Desargues transformation is in this case isogonal, [ 
and J being interchanged. Thus in the euclidean plane the absolute point 
pair should in a sense be counted with both sets of points. 

We pass now to the tetrahedron, and take as before actual perpendiculars 
on the faces as codrdinates. The equation of the absolute is then 


Q= + = 0, (7 = 1, 2, 3, 4, andz + 9), 


where Cj. is the cosine of the interior angle between the planes 1 and 2, 
and the c’s are all independent for the most general case in non-euclidean 
space. The definition for the analogue of the Fermat points is suggested 
by the dynamical conception outlined above. If four equal forces from x 
to the four vertices are in equilibrium, the bisector of the angle between 
any pair of forces bisects the angle between the remaining pair, and the 
two pairs are equally inclined to the bisector. This bisector meets a pair 
of opposite edges, and the three bisectors are mutually perpendicular. 
The configuration of four directions may be realized on a sphere by taking 
a point P and its three reflections in the vertices of a tri-rectangular spherical 
triangle. If we refer the tetrahedron to these three lines as rectangular 
cartesian axes it takes the form 


ma, 7b, 1e; pa, — pb, — pc; — oa, ob, — oc; — ty, — Tb, re. 


We are thus led to seek, as the analogue of the Fermat points, those points 
S such that the lines from S to the pairs of edges are mutually perpen- 
dicular. In euclidean space, for which our dynamical: conceptions hold, 
these points are also defined by the condition that r; + re + r3 + 14 shall 
have stationary values. The more restricted problem of making r; + 12 
+r3-+ 7, a minimum has been considered by Sturm,* who notes the 
orthogonality of the lines from such a point to the edges, but who obtains 
an analytic solution only as the intersection of three surfaces of the 12th 
order. Two problems in Wolstenholme’s collection may also be noted. 
The first of these (No. 2024) deals essentially with the dynamical aspect, 
and the second (No. 2030) with finding points such that the lines to the 
edges are orthogonal. There is no indication of any connection between 
the two problems. 


* Crelle J., 97, p. 49 (1884). 
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The three planes 


1/21, 1/22, — 1/23, —= 1/24; 1/21, = 1/22, 1/23, = 
1/21, = 1/20, = 1/23, 1/24, 


each contain two lines on the point 2 meeting pairs of edges. These 
planes are mutually perpendicular if 


2¢12 1 1 2c 
etc., (3) 

equations which are given by Wolstenholme. If now we write y = 1/z 
we get the net of quadrics 

+ yo? + = ys + ys + 2csaysys, etc., (4) 
defining a set of eight points M* which have the property that the feet of 
the perpendiculars on the faces form an equifacial tetrahedron. They are 
thus similar to the hessian points of the triangle, so that the two sets 8 


and M form a quite perfect analogue of the two-dimensional case. 
The equations of the net (4) may be rewritten 


where 7; = Cayi + + Cisys + = 1. 
The points S thus go into the points M by the collineation 
11/21 = = = 14/24. (5) 
The two sets are thus projectively equivalent. The fixed points of the 
collineation are apolar to (Q) and to the quadric 
+ + = 0 (6) 


which touches the faces at the feet of the altitudes. So far no condition 
has been placed on the c’s. If the space is euclidean Ag = 0. This con- 
dition is equivalent to the coexistence of the equations 


— Ay + + ¢13A3 + C14A4 = 0, 
C1241 — Ao + + = 0, (7) 
€13A1 + — Az + C34A4 = 0, 
€14A1 + + C3443 — Ag = 0, 


where the A’s are the areas of the faces. Comparing these with (5) we see 
that the point and plane A;, Ao, As, A, are fixed. The plane is however 
the plane at infinity, and the point is the “symmedian” point. The 


* The discovery of these is due to Dr. F. D. Murnaghan, to whom I owe in its essentials 
the above elegant treatment, which replaces my own cumbrous solution. 
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collineation is thus affine, being a pure strain whose axes are the axes of the 
quadric (6). 

The configuration of the eight points has not been in general determined. 
The equations of the net appear in a canonical form involving 6 constants, 
which is also the number of absolute invariants of a tetrahedron and quadric. 
The plane quartic obtained by equating to zero the discriminant of the 
net also appears in a canonical form, the terms in 2*y, etc., being absent. 
This would seem to indicate that the net is unrestricted, except when 
Ag = 0. The behavior of the net in special cases throws some doubt on 
this point, but the problem involves the little known subject of combinants 
of a net of quadrics, and must be postponed for the present. 

Special cases will be considered only briefly. The most important of 
these is the equifacial tetrahedron (not necessarily euclidean), defined by 
C14 = C23 = C1, ete. Writing 


= — art ae t+ as? — — + ete.; = V1 — etc., 


the quadrics 
= + = + ¢3/83 


break up into pairs of planes. Four of the points S and M coincide with the 
incenters (1, 1, 1, 1) (1, 1, — 1, — 1) GQ, —1, 1, — 1) (1, -—1, -—1,) 
and the other four form a tetrahedron in fourfold perspective to this. 
Further details of this case are omitted. It should be noted however that 
the equifacial tetrahedron, with 3 absolute invariants, gives rise to a con- 
figuration which has none; indicating the possibility that the configuration 
in the general case may not be that of a general set of “associated” points. 
Another case that can be completely solved is the “isosceles’’ tetrahedron 


C12 = Cig = Co3 = C, C14 = Co = CH = CC’. 


The net in this case contains three pairs of planes which all belong to the 
same pencil of cones, which thus have four common generators. 


UNIVERSITY OF VIRGINIA, 
March, 1921. 
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THE CAUCHY-LIPSCHITZ METHOD FOR INFINITE SYSTEMS OF 
DIFFERENTIAL EQUATIONS.* 


By L. Hart. 


1. Introduction.—The present paper is concerned with the solution of a 
denumerably infinite system of ordinary differential equations, in which the 
variables and functions assume real values, 


(1) = = [2 = 2. (21, -)]. 


Hypotheses will be imposed on (1) under which it will be proved that a 
least one solution exists satisfying given initial conditions. 

In a previous paperf the author stated certain theorems concerning 
completely continuous functions of infinitely many variables, and considered 
a system (1) in which the f; were of this type. The unique existence of a 
solution of (1), satisfying given initial conditions, was established by a 
generalization of the Picard approximation process. It was found neces- 
sary to assume that the f; satisfied a condition with respect to £ analogous 
to the Lipschitz condition used in the consideration of finite systems of 
differential equations. In the discussion which follows, the system (1) 
will be considered without the assumption of the Lipschitz condition and, 
in analogy with the corresponding result for finite systems,{ it will be 
established that at least one solution exists satisfying given initial conditions. 
The proof will be made by use of an extension to infinite systems of the 
notion of Cauchy polygons, and will be similar in method to that used by 
Montel§ in establishing the corresponding result for finite systems. 

The general system of notation used below is vectorial in character. 
After any notation has been introduced, it will retain the same meaning as 
a formal algebraic expression whenever used in the future. 

2. Extension of a Theorem of Arzela.—Let R denote the region of 
points € = (21, t2, ---), in space of infinitely many dimensions, satisfying 


the inequalities 
la; —a;| Sri = 1, 2, 


It is- well known that an analogue of the Weierstrass condensation theorem 


* Presented under a different title to the American Mathematical Society, December 
28, 1916. 

+ Transactions of the American Mathematical Society, Vol. 18 (1917), p. 125. Referred 
to in the future as Paper 1. 

t Cf. P. Montel, Annales de l’ Ecole Normale Supérieure, 3d series, Vol. 24 (1907), p. 265. 

§ Loc. cit., p. 264. 
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holds in R. That is, if S = (&,; n = 1, 2, --+) is a sequence of points in 
R, there exists a point = (21, 22, ---) and a sub-sequence (£’,) of S such 
that 

lim 2’ in = 2: [¢ = 1, 2, = Zany 

For a sequence of functions 7 = [y,(t); n = 1, 2, ---_], the notion of 
equal continuity has been defined as follows:* 

DEFINITION 1. The functions of T are equally continuous at a point to 
if, for every e > 0 a number d > 0 can be assigned so that on the interval 
|t — to| = d the oscillation of every function of T is at most e. 

If 7 satisfies Definition 1 at all points t on a closed interval (a, b), it 
is easily established that the functions of T are equally uniformly continuous 
on (a, b). In the theorem which follows there is stated, without proof, a 
property of equally continuous functions recognized originally by Arzela.t 

THEOREM I. Let the functions of T be defined and equally continuous at 
every point of an interval P which may be finite or infinite, open or closed. If 
the maxima of the absolute values of the functions yn(t), for t on P, possess a 
common finite bound, then, we may extract a sub-sequence from T which con- 
verges at every point of P. Moreover, the convergence is uniform on every 
closed sub-interval of P. 

Consider a sequence S’ = [£’,,(t); n = 1, 2, ---_] where, for every ¢ on 
P, é’,(t) isa point in R. Let the sequence formed by the 7th codrdinates of 
the functions of S’ be represented by S; = [2’in(t); n = 1, 2, --- ]. 

TueorEM II. For every i suppose that the functions of S; are equally 

continuous on P. Then, we may select a sub-sequence S = [En(t) | of 8S’ and 
a point y(t) belonging to R for every t, such that, for ail t on P, 
(2) lim vin(t) = (i = 1,2, ---), 
where tin(t) and y,(t) are the i-th codrdinates of —n(t) and n(t) respectively. 
Moreover, the convergence in (2) is uniform on every closed sub-interval (a, b) 
of P. 

On applying Theorem I to the sequence Sj, it is seen that we may select 
a sub-sequence S® = [£,“(¢) ] from S’ corresponding to which a function 
yi(t) exists satisfying 
(3) lim = 


uniformly on (a, b), where 21,)(¢) is the first coérdinate of —,(#). As 


_aresult of applying Theorem I to the sequence formed by the second 


coérdinates of the functions of S®, it follows that a sub-sequence 


*Cf. Montel, loc. cit., p. 236. 
+ Cf. Montel, loc. cit., p. 237. 
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S® = [é,(@] of S® and a function y2(f) may be selected for which 
there holds 
(4) lim ain (t) = yill) = 2, k = 2), 


uniformly for ¢ on (a, b). Obviously, because of (3), condition (4) is also 
true when 2 = 1. Ina similar manner, for every h we obtain a sequence 
= and a function y,(t) such that the codrdinates of 
satisfy (4) with (k = h; i= 1, 2, ---, h). It is easily established that a 
sequence § satisfying Theorem II is obtained by placing &,(¢) = &(d), 
and by identifying y;(¢) of (2) with that in (4). In view of (2), it is clear 
that every y;(t) satisfies |y;(t) — a;| = ri. 

3. The Existence Theorem for System (1).—In equations (1) suppose 
that (é, ¢) isin the space U defined by the condition ( in R, a St SD). 

DEFINITION* 2. A function f(£), defined in R, is completely continuous 


at a point in R af, whenever limn=» Vin = 1, 2, -+-), follows that 
lim f(&n) = fé [En — (Xin; Lon, °° -), in R}. 


It is well known that a function f(), completely continuous at all points 
of R, possesses a maximum and a minimum, each of which is attained at 
least once in R. If, for every t, &(t) is a point in R whose codrdinates 2;(t) 
are continuous functions of ¢, it is easily verified that f[£(¢)] is a continuous 
function of t. The author has provedt the simple result that, if (51, be, ---) 
is a sequence of positive numbers and if e > 0 is assigned, a number d > 0 
can be found such that if (¢’, £’”) are two points of R whose ith codrdinates 
satisfy |a;’ — x/’| = db; (i = 1, 2, ---), then 


(5) — FE) | Se. 


In the future let us assume that f; in (1) satisfies the following conditions: 
Hy. If [(én, tn); nm = 1, 2, «++ ] is a sequence of points in U for which 
lim = lim 23n = 2; (2 = 1, 2, --+), 

it follows that lim,—, fi(En, tn) = fi(€, 2). 

H,. There exists a number M > 0 such that, for every value of 7, the 
maximum m; of |f;(£, ¢)| in the region U satisfies m; = r;M. 

THEOREM III. Suppose that aSto<b. Then, there exists at least 
one function &(t) whose codrdinates x;(t) fulfill the initial conditions (x;(to) = ai} 
a= 1,2,---) and satisfy (1) for t) = t S c where c is algebraically the smaller 
of b and (tp) + 1/M). 

In order to prove the theorem let us extend the notion of a Cauchy 


* Cf. Paper I, p. 129. 
Paper I, p. i380. 
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polygon,* as defined for a finite system of differential equations, to the case 
of an infinite system (1). 

Let p represent a partition of the values ¢ = ¢) with the division points 
(to, tm = 6). Then, the codrdinates 6,(t) of the Cauchy polygon 
B(t) for (1), corresponding to the partition p and satisfying B(to) = a 
= (a, dz, +++), is defined by the following equations, where i = 1, 2, - 


b(t) = + fila, to) (¢ to) (to = t = th), 
b(t) = + [hStSh; B= 


(6) a 
b(t) = bith) + fi(Bas — te) [th St S tars; Ba = Br) J, 


In order to obtain the sequence of polygons which will later be shown to 
converge to a solution of (1), select a sequence of partitions (p’n; n = 1, 2, 
.--) with norms (d’,). Suppose that lim,—. d’, = 0. To each partition 
p', there corresponds by (6) a polygon §’,(¢). Theorem III will be proved by 
showing that the sequence S’ = (8’,(¢)) satisfies the hypotheses of Theorem 
II and by then establishing the fact that the limit function given by that 
theorem is a solution of (1). 

First consider some properties possessed in common by all Cauchy 
polygons. For every partition p, it follows from (6) and Hz that, if ¢ is 
on the interval (¢;, t,41), where = ¢, then, 


|b:(t)—a;| = | b:(t) — bi(th) |b — (ths) i+ | bi(t1) — bi (to) | 


Consequently, for all ¢ on (to, c), the polygon B(f) is defined and in R. In 
the future suppose that ¢ is on the interval (fo, c). In the same manner as 
we obtained (7) it is verified that, for every partition p, and for all points 
(t’, t’’) on interval (fo, c), 
(8) — b:(t”’)| S Mr; — t’’|. 
Hence, if e > 0 is assigned, a number w > 0 can be found so that, for every 
partition p and for every 2, 
(9) |b,(t’) — b(t’) | Seri, 
provided that |t/ — ¢’’| w. 

Consider the sequence 8S’. As a consequence of (7) and (9) it follows 


* 7. Bliss, Princeton Colloquium, p. 89. 
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that S’ satisfies the hypotheses of Theorem II on the interval (fo, c). Let 
S = (6,(t)) be the sub-sequence of S’ and let &(¢t) be the point given by 
Theorem II in the present instance. Equation (2) becomes 


lim Din(t) = LBn(t) = (din(t), bon(t), 7= 1, 2, 


uniformly for ¢ on (to,c). Asa consequence of (6) it is obvious that £(t)) =a, 
In view of (9) it is seen that, if e > 0 is assigned, a number w > 0 may be 
found such that, if |t/ — t’’| S w, then 


(10) | a(t’) a,(t’’) | = Or; = 2; 
It remains for us to show that &(#) is a solution of (1). 

As a consequence of (6) we may state that the codrdinates b;n(t) of 
B,(t) satisfy 

bin(t) = hy + Bnd), t jdt, 

where ¢ is a function of the variable ¢ of integration and is defined as the 
last division point preceding ¢ in the partition p, corresponding to 8,(t). 
Let us show that, for every value of 7, 


(11) lim f[8.(@, #] = 
uniformly for ¢ on (¢, c). When this has been proved, we may write 
a(t) = ait SiflEO, tldt, 


from which it follows by differentiation that &(¢) is a solution of (1). 
To establish (11) let us consider 


(12) = t | t || + t | t’]| 
+ lf Brit’), t’ fLEC), t’ + fe’), t’ fi E(é), t ||, 


where ¢’ is any point on (to, c). In the function f;(x1, 22, ---; ¢) we may 
think of (a1, a, ---; ¢) as being the coérdinates in the space of infinitely 
many dimensions defined by (|a; — Sri; to Stc). As a conse- 
quence of (5), if e > 0 is assigned, a number w > 0 can be determined so 
that, if — and |t’ Sw, 


In (12) recall that |¢ — ¢| S d,, the norm of the partition p,, and, therefore, 
as a result of (8), , 

| bin(t) bin(t) | d,Mr;. 
Therefore, it follows from (13) that the first term on the right in (12) 
approaches zero for n = ©, uniformly for all ¢ on (to, c). In view of (9) 


J 
a 
| 
a 


Hart: Infinite Systems of Differential Equations. 231 


and (10), an application of (13) to the second and fourth terms on the 
right of (12) shows that, if ane > 0 is assigned, a number g > 0 can be found 
so that, if {tf — ¢’| S g, then each of these terms will be at most e/3, for all 
values of n. 

Before considering the third term in (12), let us form a partition (v» = fo, 
V1, «+, Um = ¢), of the interval (fo, c), with the number g of the previous 
paragraph as its norm. Since f;(&, ¢) is completely continuous and since 

lim b;,(t) = 2;,(t), 


an integer VN may be determined so that 


If LBn(v;), 05] — 9,]| Ss G=0,1,---,m; 
In (12), for a fixed value of ¢t, suppose that t’ is one of the (v;) differing 
from t by at most g. Then, if n = N, the sum of the last three terms in 
(12) is at most e. Since N was chosen independently of ¢, it is seen that 
the uniformity stated in (11) has been completely established. Hence 
Theorem III has been proved. | 

Theorem III was concerned with the existence of a solution on the 
interval t 2 A similar thoerem may be stated for the interval a St S fy. 

In paper I the author established* the existence of a unique solution of 
(1) under the conditions of Theorem III together with an additional assump- 
tion. The proof of the uniqueness could be made in the same way in the 
present paper with the results of Theorem III as a starting point. 


University OF MINNESOTA, 
July 1, 1920. 
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BOUNDARY VALUE AND EXPANSION PROBLEMS; FORMULA- 
TION OF VARIOUS TRANSCENDENTAL PROBLEMS.* 


By R. D. CARMICHAEL. 


Proceeding under the guidance afforded by a previous investigation 
of the algebraic basis of the theory of certain boundary value and expansion 
problems we here develop the formal aspects of several transcendental 
expansion problems which are rather different externally but are abstractly 


of much the same essence. In each case the work is carried forward to the 
point of derivation of the form of the expansion and the values of the 
coefficients; and the problem of determining the nature of the functions 
having expansions of each of the kinds enumerated is thus clearly set forth. 
It seemed desirable to present in a single paper this formulation of several 
problems, so that the reader may most conveniently compare their elements 
of similarity or difference; and to leave to further investigation the detailed 
solution of each problem deserving such further treatment. It will be 
apparent to every one that we have given merely a selection of problems 
from a much larger number which emerge in a similar way from the con- 
sideration of various sorts of limiting cases of the algebraic systems dealt 
with in the paper already mentioned. 

In § 1 we treat the boundary value and expansion problem arising in 
connection with a single linear differential equation involving a single 
parameter. In § 2 the corresponding problem for a system of first order 
equations is more briefly treated and the results are carried over by 
Volterra’s limiting process so as to yield the formulation of a similar expan- 
sion problem for certain integro-differential systems; and in § 3 the results 
of § 2 are rapidly generalized to problems involving vy parameters. Various 
problems involving difference and integro-difference equations are treated 
in § 4; these are analogous to the corresponding ones for differential and 
integro-differential equations developed in §§ 1, 2, 3. In §5 is given a 
brief statement of a great variety of mixed problems involving one or more 
parameters. The character of the problems arising from partial differential 
equations and integro-differential equations is developed rapidly in § 6 
and certain typical cases are briefly treated separately. A particular case 
of problems arising from partial differential equations is one which is con- 
nected with the theory of vibrating plates; this is briefly treated in § 7. 
 * Presented to the American Mathematical Society, April, 1920. 
+t AMERICAN JOURNAL OF Matuematics, Vol. XLIII (1921), p. 69. 
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More details are given in § 1 than elsewhere in the paper, the justification 
being that those which are omitted are fairly obvious to one who has in 
mind the corresponding developments in § 1. 

1. The Problem Arising from a Single Differential Equation with One 
Parameter.—In connection with the linear differential expressions 


d™u 


m—1 
dx dz” dz 


let us consider also the Lagrange = expressions 


1” 
M(v) = (- [lnv] + (— — Che] + lov 
= My ~+ Mov, 
M,(v) = ., ~+ hov. 
dx™ 


Here it is understood that m may be zero so that L,(w) may have the par- 
ticular value gow. We assume that the coefficients J;, gi, mx, hi, with 
their derivatives of all orders, are real-valued single-valued and continuous 
functions of x in an interval (ab) and that J, and gm do not vanish in (ab). 
Without loss of generality we take /,, to be positive in (ab). With the linear 
differential equation 


(1) L(u) + AL, (u) = 0 

and certain n linear homogeneous conditions in u(a), w’(a), «++, (a), 
u(b), u’(b), 

(2) Ui(u)=0, = 0, 

we associate the adjoint linear differential equation 

(3) M(v) + AMi(v) = 0 

and n like adjoint conditions 

(4) Vasi(v) = 0, V n+2(v) = 0, Von(v) = 0 


in such wise as to generalize the corresponding problem treated by Birkhoff* 


for the case when L,(u) = gow. 
We must first make precise the boundary conditions. In doing this 
as well as in developing the other parts of this section we follow the sugges- 


* Transactions of the American Mathematical Society, 9 (1908), 373-395; 219-231. 
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tions afforded by Birkhoff’s treatment of the special case. If we write 


n—1 


d 


d d? d™y 
| Ino — © (lav) + | 


ly 


d 
y)* = eee — —— (Gor) y) —— 
(— ie (Gm) 
then it is well-known and easily verified that we have the identities 
(5) vL(u) — uM (rv) = {R(u, v)}, vL;(u) — uM,(v) = : {Ry(u, v)}. 
x dx 
Let us now consider the quantity 
[R(u, v) + ARi(u, v) 


as a bilinear form in the two sets of variables w‘"—) (a), uw" (a), «++, u’(a), 
u(a), (b), (b), u’(b), u(b) and v’(a), ---, (a), v(b), 
v'(b), «++, oY (6), and let the form be arranged in a square array whose 
columns contain the w’s in the order written and whose rows contain the 
v’s in the order written. The matrix of this bilinear form obviously has 
the following properties: every element below the main diagonal is zero; 
every element in the upper right hand fourth of the matrix is zero, the 
division being made by horizontal and vertical lines through the center; 
the determinant of the matrix has the non-zero value {/,(a)/,(b)}". 
Now this bilinear expression may be written in the normal form 


2n 


(6) [R(u, v) + AR(u, v) = U;(u)V ;(v), 


where the U; and the V; are linear and homogeneous functions with con- 
stant coefficients, the first in the variables uw and the second in the variables 
v. From the totality of such representations we propose to choose those 
having certain properties. 

In case L,(u) = gow, so that our bilinear form is independent of \, 
the linear forms U; may be chosen in any way so that the 2n forms are 
linearly independent. The forms V’; are then uniquely determined and are 
linearly independent. This is the choice to be taken in the case of the 
Birkhoff problem. 


*Tf m = 0 we take R,(u, v) to be identically zero. 


& 
: 
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In the more general case when the bilinear form contains \ we wish to 
realize the following conditions: the linear forms U;, Us, ---, Un, Vn, 
Von Shall be independent of A; the whole set Ui, Us, ---, Usn 
shall be linearly independent for every value of \; the whole set V1, Vo, 
.++, Von shall be linearly independent for every value of X. For the last 
two conditions it is necessary and sufficient that the determinant of the 
set of linear forms U; shall be independent of d (except possibly for a factor 
dependent on ) but vanishing for no value of \) and different from zero, as 
one sees readily by aid of the fact that the determinant of the U; multiplied 
by the determinant of the V’; has a value equal to the (non-zero) value of the 
determinant of the bilinear form. That these three conditions may in 
fact be realized may be seen from the following particular choice of linear 
forms: 

U,=u* U.= ua), ---, wa), = u? 6); 

= (a), = 0'(b), = 0''(b), +++, Von = (6); 


for 2 =n, V; = total factor of U; in the fundamental bilinear form; for 
i>n, U; = total factor of V; in the fundamental bilinear form. Indeed 
it is now evident that the needful properties may be realized in a great 
variety of ways. 

When the forms U; and V; in (6) are determined in accordance with 
the conditions named, then equation (1) and conditions (2) set a given 
boundary value problem the adjoint of which is set by equation (3) and 
conditions (4). It is clear, moreover, that problem (1), (2) is likewise the 
adjoint of problem (3), (4), so that the relation between the two problems is 
a reciprocal one. 

Remark 1.—We have seen that the conditions named can always be 
realized whatever the value of m, subject merely to the condition that it 
shall be less than n. For smaller values of m we also have certain addi- 
tional freedom. We may more generally take & linear forms U; and 2n — k 
linear forms V; independent of \ where / is not necessarily equal to n but 
may be less or greater by an amount depending on m. When m is zero k 
may be any positive integer less than 2n.* We may then make correspond- 
ing modifications in the two adjoint boundary problems. For the purposes 
of this paper, however, it is desirable to realize the maximum symmetry in 
the relations of the two problems; and consequently we confine our atten- 
tion to the case indicated in the definition as explicitly stated above. 

Remark 2.—It is proper to inquire whether the restriction that m shall 
be less than n is essential. The case m > n reduces to the case m < n 
by replacing \ by 1/A (except for the particular value \ = 0), so that we 
 *See a treatment of this case by Bécher, Transactions of the American Mathematical 
Society, 14 (1913), 403-420. 
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need to consider only the case m = n. In this case the determinant of the 
bilinear form in the first member of (6) is not in general independent of ); 
and hence the properties of the forms U; and V; could not be maintained in 
general, properties which are needed in the further development. These 
properties might, however, be retained in large measure when m= n 
provided that g, and its derivatives have suitable values at a and at b; 
but we shall not attempt to carry this case along with the other, owing to 
the special hypotheses which would be necessary and the complications 
which would raise in certain parts of the work. 

Remark 3.—The question of generalization to an equation of the form 


L(u) + ALy(u) + + +++ + ML, (u) = 0 


also arises. Under suitable hypotheses it is not difficult to maintain the 
necessary characteristics of the boundary conditions. Some of the theorems 
are still maintained without change; but as early as theorem III certain 
difficulties would arise in general with this more extended case. 

We may now readily prove the following theorem: 

TuEorEM I. Jf for \ = a solution* ii(x) of (1), (2) exists, then a solu- 
tion 0(x) of (3), (4) also exists for \ = d; uf %(x) is unique (except for a con- 
stant factor), 0(x) 7s also unique (eacept for a constant factor). 

Since @ is a solution of (1), (2) for \ = A, we have 


U,(a) = = --- = U,@ = 0, 


and for some k 
0, 


where in U,4; we have the value X of X, since otherwise we should have 


= u' (a) = (a) = = ii’ (b) = eee = (b) = 0 
and our function (2) would be identically equal to zero, contrary to hy- 
pothesis. 


In the n-fold linear spread of solutions v(x) of (3) for \ = X, there is at 
least one, say o(”), which satisfies the n — 1 linear homogeneous conditions 
= 0, 
where 2 runs over all the numbers of the set 1, 2, ---, m except k. Now 
we have L(%) + ALi = 0 and M(o) +AM,(0) = 0. Hence if we mul- 
tiply the second identity in (5) by \ and add to the first we find on in- 
tegrating that + 3) is a constant. Therefore from (6) we 
see that Unix (%)Vnir(%) = 0; whence it follows that V,,,(%) = 0, and 

* By “solution” we usually mean a function which satisfies the equation and condi- 


tions but is not identically equal to zero; when we intend otherwise we shall indicate that 
fact explicitly. 
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hence that all the conditions (4) are satisfied by 3, so that it is in fact a 
solution of (3), (4). 

It remains to be shown that 7 is unique whenever &% is unique. For 
this purpose let us suppose, if possible, that is not unique. Then two 
linearly independent solutions of (3), (4) exist, say 7 and v*. Then k and 7 
exist such that 


Vian, 
for, otherwise, constants c and d would exist,* not both zero, such that 
V (cd + dv*) = 0 for 1= 1, 2, ---, n, whence it would follow that the 
solution v = cd + dv* would satisfy all the conditions V;(v) = 0, 2 = 1, 2, 
..», 2n, and this is impossible since the V’s are linearly independent for 
every value of \. Now choose w*, linearly independent of i, so as to satisfy 
the n — 2 conditions 


U;(u*) = 0 


for i running over the numbers of the set 1, 2, ---, , except k andj. Then 
from (5) and (6) we see that we have 


+ U;(u*) V;() = Q, 
+ Uj(u*)V;(0*) = 0. 


From this we see that U;(u*) = 0 = U;(u*). Hence u* satisfies (1), (2). 
Therefore, if the solution of (3), (4) is not unique for \ = d neither is the 
solution of (1), (2) unique. Hence we conclude to the truth of the second 
part of the theorem. 

A value of \ for which the system (1), (2) [and hence the system (3), 
(4) ] has a solution will be called a characteristic value. The characteristic 
value is said to be simple if the solution corresponding to it is unique (except 
for a constant factor). 

THEOREM II. If yi, yo, +++, yn are n linearly independent solutions of 


(1) for \ = X, the condition that is a characteristic value is that the deter- 


*In more detail one may proceed thus: distinct numbers k; and j; exist such that 


Vn@) #0, +0. 
If then — 

Vn = Vino) Vi) 
we have 


Vio") $0, Vigo) $0. 
But if our determinant is always zero it is sufficient to take 
c = V;.(v*), = — Vi,(r), 


in order to realize the result stated. 
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minant 
+++. Uilyn) | 
Us(y1) Us(y2) +++ 
if (y1) U, (y2) (Yn) 


shall vanish; a characteristic value d of d is simple when and only when some 
jirst minor of A does not vanish. 
If we take the general solution of (1) in the form 


U = C1y1 C2Y2 treet CnY ns 


we see that the vanishing of A is a sufficient condition that conditions (2) 
may be satisfied through an appropriate choice of ¢1, ¢2, -*+, ¢n. When 
some first minor does not vanish this choice is unique (except for a constant 
factor). 

THEOREM IIT. Jf u;(x) and v;(x) are solutions of (1), (2) and (8), (4), 
respectively, the first for \ = d; and the second for \ = dj, and af dx + dj, 


then we have 
= 0, = 0. 


From the equations 


L(u;) + ALi (ui) = 0, M(v;) + = 0 
we have 
L(ui) + ALi (us) ] — ud M(o;) + J} + Oi — Au:Mi(o;) = 0. 
If we integrate from a to b we see through aid of (2), (4), (6) that the expres- 
sion inclosed in braces yields the value zero. Hence we have 


(Ai — Aj) = 0. 


Dividing through by the factor \; — ; we have the first relation in the 
theorem. The second may be proved in a similar way. Moreover, either 
of these relations follows readily from the other by means of the second 
relation in (5) and relation (6), looked upon as an identity in A, together 
with the boundary conditions. 

If we have two sets of functions w(x), W(x), u3(v), and v(x), 
v3(2), and the associated adjoint linear forms L,(u) and M,(v) such 
that relations (7) are satisfied for every i and 7 which are different we shall 
say that the two sets of functions are adjoint each to the other. If the two 
sets are the same and the form L,(uw) is self-adjoint we shall say that the 
single set is self-adjoint. The conceptions thus introduced are obviously 
generalizations of the conceptions of biorthogonal sets and orthogonal 
sets of functions. 


If now we have further 


| 

| 

| 


it is clear that we may formally determine the coefficients c; in the formal 
expansion of a given function f(x) in the form 


f(x) = (x) + + 


In fact, if one multiplies by J/,(v;) and integrates from a to 6 one has for- 
mally 


When J is not a characteristic number for (1), (2) there exists* a unique 
G(x, s; 4) such that the solution ¢ of 


L(y) +ALi(y) = 17, = U2(y) = = Unley) = 9, 


is given by 


g= fG(a, 8; r)r(s)ds; 
and likewise there exists a unique H(x, s; \) such that the solution y of 


Mi) = = Vasey) = = Von(y) = 0, 
is given by 
Y= (a, d)r(s)ds. 
Moreover, we have 
G(x, 8; X) = H(s, 2; d). 
Explicit formulae for G and H are readily obtained (see Bocher, l.c.). 
The formula for G may be written in terms of a fundamental system 4, ye. 
-+, Yn of solutions of (1) as follows: 


yi(z) +++ 8; d) 


N(a, A) _ U s(n) Unlys) U, un) 


(8) G(a,s; A) = A(A) 


where A(A) has the value given in theorem II and G(z, s; \) has the value 


Y2(x) Yn(2) 


1 Un (8) 


2 In (8) yo™D(s) | 

the sign + being taken when x > s and the sign — when z < s. 
* Boécher, “Les Méthodes de Sturm,”’ Chapter V. 
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Now the functions y, y2, +++, yn may be taken analytic in X. Hence 
G is analytic in \ except when A(A) = 0, that is, except for the characteristic 
values \; of X for problem (1), (2). If A; is a simple characteristic value 
we may write 
s) 


(10) Gla, 8; = 0) 


where s; A) is analytic in \ = d;. From (8) we have 
fm a) 83 
R;(2, S) = 


where A’(\;) denotes the value at \ = \; of the derivative with respect tod 
of A(A). In N(a, 8; d,) the coefficient of G is zero. Hence N(z, s; },) 
and its first n derivatives with respect to 2 are continuous. Moreover, it is 
easy to see that it satisfies (1), (2); therefore R,(x, s) is a solution of (1), 
(2). From the relation between G and H we infer that R;(x, s) is a solution 
of (3), (4) ins. From the definition of simple characteristic value it follows 


then that 
R;(x, 8) = cu;(x)0;(s) 


where c; is independent of x and s and is yet to be determined. 
Now from (10) we have 
lim {(A — A,)G(a, 8; A) — eui(x)v;(s)} = 0; 


A=Ai 


whence it follows that 
(11) 8; — exui(x) Li (ui(s))ds_] =0. 
A=Aj 


But since 
L(ui) + (ui) = A Ui(u;) = U2(u;) = +++ = Un(ui) =0, 


we have from the fundamental property of G that 
(A — Ai) (a, 8; A)Li(ui(s))ds = 
Substituting this value into equation (11) we have 
Sa'vi(s) Ly (ui(s))ds = 1. 
Similarly we should obtain the relation 
(s) My (0;(s))ds = 1. 
Hence we have for the residue R;(z, s) the value 


ui(x)v;(s) _ ui(r)oi(s) 


If I is a contour in the \-plane which incloses (just) the characteristic 


(12) R,(x, 8) = 


ix 


nee 
Stic 
ue 
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numbers Ai, Az, -**, An, and if the latter are simple, we have readily 


M; 
= 


f(s)dsdr 


V— 
1 My (v;(s))ds 


Thus we have represented in the form of a contour integral the sum of 
a finite number of terms of the formal expansion 


Sa’ f(s) Mi(o,(s))ds 

(14) f(x) = (8) (s)) U;(2) 
of a function f(a) for the case when the characteristic numbers Aj, 2, Az, 
- are all simple and G(a, s; \) has a pole of the first order. At a value 
\; of \ for which these conditions are not satisfied, the corresponding term 

in the formal expansion in (14) should be replaced by 


(ce, 8) 


where R;(z, s) denotes the residue of G(x, s; \) at X = X; and 2;(s) is the 
solution V(a, s; ;) of (3), (4) ins. Thus we have a definite formal expan- 
sion for every given function f(a) for which the integrals involved exist. 

In accordance with the general plan of this paper we rest this investiga- 
tion at this point. We have carried it far enough to set clearly the funda- 
mental expansion problem associated with it, namely, the problem of deter- 
mining those functions f(a) for which the given formal expansions are valid. 
This problem Birkhoff (l.c.) has treated with great success for the special 
case considered in his paper. It appears likely that the methods of 
Birkhoff’s paper are sufficient to carry through the more general investiga- 
tion. The results obtained above are sufficient to initiate the further 
treatment. 

In the following sections we derive formal expansions in connection with 
a great variety of boundary value and expansion problems. We do not 
go into as much detail as in this section; but the reader who has this section 
in mind will have no difficulty in seeing how the results actually obtained 
in each case set the fundamental expansion problem for that case. 

2. The Problem Arising from a System of Differential Equations of the 
First Order with One Parameter and from Related Integro-differential 
Equations.—Let us consider the adjoint systems of differential equations 


dy; 
dx 


j= 


Qj; 


| 

is 
1), 
on 
Ws 
0. 
0, 

@) dx ( 
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If for fixed 7 we multiply these respectively by z; and y;, add the resulting 
equations member by member, sum as to 7 from 1 to n, and then integrate 
from a to b, we have 


Let us now suppose that homogeneous linear boundary conditions 
(involving the end-points a and b) are set up on the y’s and adjoint condi- 
tions on the z2’s so that (3), for the solutions to be considered, is satisfied in 
virtue of the boundary conditions alone. Moreover, let us suppose that this 
is done in such a way that the characteristic values of \ for (1) with its 
conditions are the same as those for (2) with its conditions, the situation 
being analogous to that treated in § 1, and that the number of characteristic 
values is infinite. Let us denote the distinct characteristic values and 
corresponding solutions by the symbols \™, y;, 2;, k = 1, 2, 3, ++: 
In (1) let us now replace \ and y; by \“ and ats in (2 2) let us replace \ 
and 2; by A and 2z;, where \ +)”; then we have two new systems. 
For fixed i let us multiply the former by 2, and the latter by y;“, add, 
sum as toz from 1 to n, integrate from a to b, and simplify the first member 
by means of the equation to which (3) reduces when y; and 2; are replaced 
by ys and 2;. Thus we have 


= SX, + AM axis) yj — $1; 
whence, on simplifying and dividing by A — \, we have 


(4) = 0, 


f= 


If we suppose that the integral in the first member of (4) is different in 
value from zero when k = 1, we have a convenient means for the formal 
determination of the coefficients c, (independent of x and ») in the simul- 


taneous formal expansions of n given functions f(a), fo(~), fn(x): 
&=l 


For determining these constants c, we proceed as follows: multiply both 
members of (5) by a;:z;(, sum as to 7 and as to j from 1 to n, and integrate 
from a to b; in this way we obtain a relation from which we have readily 
the value of c;, namely: 


dx 
(6) = = imi. ..., 
Ode 


i=1 j= 
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Thus we have in (6) the formule for the coefficients ¢;, the same for each 
function f;(x), needed in the joint expansions in (5) of the n functions f,(x). 
They are formule of great interest. 

We have outlined very briefly the foregoing expansion problem, which 
is probably the same as that of an unpublished paper of Birkhoff (see 
abstract in Bulletin of the American Mathematical Society, (2) 25 (1919), 
p. 442), because we wish to utilize the suggestion afforded by it for setting 
up a related problem for certain integro-differential systems obtained from 
the foregoing differential systems by the limiting process which Volterra 
has so frequently employed and to such great advantage. It is clear that 
the notion of joint expansion utilized in connection with equation (5) is an 


_ essential element in the formulation of the problem. This notion I obtained 


not from Birkhoff’s abstract (which is too brief for details of any sort) but 
from the dissertation of Dr. C. C. Camp which the author was kind enough 
to allow me to read in manuscript.* In this dissertation there is a detailed 
and full treatment of a very special case of the expansion problem formulated 
above; and the range of ideas developed has been of great use to me in 
the formulation of the boundary value and expansion problems of this paper. 

If we apply to systems (1) and (2) the Volterra limiting process we may 
obtain the integro-differential equations 


(8) = {— u(a, t, — dv(a, t, 8) }2(a, 


where the range of variation of x is from a to b while that of s and ¢ is from 
ato 8. Let us multiply both members of (7) by z2(2, s) and both members 
of (8) by y(a, s), add the two resulting equations member by member, and 
integrate with respect to s from a to 8; then, on carrying out certain 
simplifications through an interchange of the order of integration with 
respect to s and ¢ and through certain obvious reductions, we have formally 
the relation 


0 
= tyla, 8) }ds = 0. 


Integrating with respect to x from a to b, we have 
(9) SF {y(b, 8)2(b, s) — y(a, s)z(a, s)}ds = 0, 


a relation which corresponds to (3) for the differential systems. 


* Since this was written a paper by A. Schur has appeared dealing with this problem; 
see Mathematische Annalen 82 (1921): 213-236. It contains (p. 216) a reference to a 
special case of the problem treated by Hilbert (Géttingen Nachrichten, 1906, pp. 474-480). 
This I had not seen before. 

Dr. Camp’s article will appear in the next volume of this JournaL. [Editor.] 
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The first member of the last equation is a sort of bilinear form in a non- ; 
enumerably infinite set of variables y(a, s), y(b, 8), z(a, 8), 2(b, s), a set of four 
variables being formed for each s of the interval (a8). One desires to have 
a reduction of this to a normal form analogous to that used in connection 
with equation (6) of §1. At present we do not attempt to solve this 
problem in general. What we desire is to have a set of boundary conditions 
on y and z separately, linear and homogeneous in character, so that these 
boundary conditions alone shall imply equation (9) and so that the y- 
problem and the z-problem shall have the same characteristic values \ 
infinite in number. Assuming that such equations have been set up, we 
seek to determine the formal character of the boundary value problem in 
its simplest aspects. 

Let us represent characteristic numbers and corresponding solutions 
by Ax, ye(a, 8), 2x(2, 8), k = 1, 2, 3, ---. Let us replace \ and y in (7) by 
A, and yz; let us replace \ and z in (8) by A; and 2), k + /; let us multiply 
the first resulting equation member by member by z;(2, s) and the second 
by yx(x, s), add the two resulting equations member, integrate with respect 
to s from a to B and with respect to x from a to b, and simplify the resulting 
relation by aid of what equation (9) becomes when y and z are replaced by 
y, and 2, respectively. From both members of the equation which so 
results remove the non-zero factor A, — 7; thus we have 


(10) See SF SF t, 8)yx(x, 8)zi(a, t)dsdtdr = 0, + 


If we suppose that the first member of (10) is different from zero for 
k = land if we seek constants c; so that a given function f(x, s) shall have 
the formal expansion 


(11) fiz, 8) = CrY K(X, 8) 


we are led (through the relation obtained on multiplying both members of 
(11) by o(a, t, s)z:(x, t) and on integrating appropriately) to the formule 
for c;, namely: 


SESE v(x, t, 8)yr(x, 8)2(x, t)dsdtdx 


As a special case of the foregoing problem let us take that in which u 
and v are independent of x and let us write 


u(x, s, t) = g(s, 3, t) = h(s, 
Solutions exist of the forms 


y(x, = p(x)o(s), 8) = R(x) S(s), 


or 
ve 


of 
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where each function in the second members is a function of the single 
variable indicated. As special boundary conditions implying equation 
(9) we may take the following: 


y(b, 8s) — y(a,s)=0, 8) — 2(a, 8)'= 0. 
These conditions now reduce to the following: 
p(b) — p(a) = 0, R(b) — R(a) = 0. 


If we substitute the foregoing special values of y and z into the special 
forms of (7) and (8) indicated we may separate variables in the way which 
is classic in the study of the partial differential equations of physics and so 
come through to the following relations: 


(13) p’(x) = up(z), p(b) — p(a) = 0, 
wa(s) = {9(s, t) + dh(s, t)}o@de; 


R’'(x) = — pR(a), R(b) — R(a) = 0, 


(14) — uS(s) = t) + AR(s, 


where yw is a constant. In (13) and (14) we use the related constants pu 
and — uw so as to maintain the requisite property of adjointness. The 
differential system in each of these pairs of conditions is easily solved, 
appropriate characteristic values of ~ being thus determined. Then the 
integral equation of the pair is to be solved, appropriate characteristic 
values of \ having been determined for each value of y. 

Let us further specialize. Suppose that we take b = 7 and a= — zr. 
Then for » we have the values n V— 1, where n is a positive or negative 
integer, and for p(x) and R(x) the corresponding values e”*"-! and e~"*'=1! 
respectively. If we take g(s, t) = 0 and employ a non-zero value of u 
we have in the second line of (13) and in the second line of (14) an integral 
equation with parameter \/y in the one case and — X/y in the other case. 
It is well-known that the quantities a, 8, h(s, t) in this equation may be 
so chosen that for each value of uw corresponding values of \ exist giving 
rise to the functions cos ms for solutions where m runs over the set 1, 2, 3, 

For the special case thus defined the expansion (11) has the form 


f(a, s) = (Cmn COS + Ymn COS mse~”"*), 
m, n==1 

If f(z, s) is a real-valued function when x and s are real it is clear that 

Cmn and Ymn Must be conjugate imaginaries. Then the expansion reduces 

to the form 


f(x, s) = 2 (Amn COS NX + bmn SiN NX) Cos ms, 
l 


n, r= 
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where @mn and bm, are real numbers. Furthermore, if f(z, s) is an even 
function of x as well as of s we must have b,,, = 0 so that the expansion 
becomes the double cosine series 
f(x, 8) = amn cos ms cos nz. 
m,n=1 

Thus we see that the expansions afforded by (11) and (12) contain as special 
cases some of the classical expansions of analysis. 

The problems formulated in connection with (1), (2) and (7), (8) are 
capable of various generalizations. It will be sufficient to indicate briefly 
one of them. Let us consider the system of equations 


(15) #) Dd (aij + Aas) 8) 
Ox Jj=1 


n 


+ 8, t) + Ao t) Jy; (x, t)dt, 


(16) = (— — ays) 8) 


+ [~ t, s) t, 8) ]z;(2, t)dt, 
j=l 


for? = 1,2, ---,m. Here in general the a;; and a;; are functions of 2 and 
s. In the special case when p;; = 0 = o;; and aj; and a;; are functions of x 
alone, we may take y; and z; to be functions of x alone, and we have our 
systems (1) and (2). It is clear that equations (7) and (8) are included in 
another way as a special case. 

Let us multiply (15) by z:(2, s) and (16) by y,(a, s), integrate the resulting 
equations as to s from a to 6 and in the relation coming from (16) interchange 
the order of double integration and in the part of the formula containing 
this double integration interchange s and t, add the two equations thus 
gotten from (15) and (16), and in the resulting equation sum as to 2 from 
1 to n. Finally in this resulting equation integrate as to x from a to b. 
Then we have the relation 


(17) SELL s)z;(b, 8) — y;(a, s)z:(a, s)} = 0. 


This corresponds to our previous relation (9). 

Let us suppose that linear homogeneous boundary conditions on the 
y; and the z; separately are so chosen that relation (17) is implied by these 
boundary conditions, that the y-problem and the z-problem have the same 
characteristic values \ infinite in number, and let us denote the correspond- 
ing values and solutions by dz, y:, 2;, k = 1, 2, 3, -- Then by a 


ven 
sion 


dt, 
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method of procedure which is now obvious we come through to the funda- 
mental reciprocal relation 


s)yi (a, 8)z; (a, 8) ]ds 
(18) 
+ fF t, (x, 8)2z; (a, t)]dsdt}dx = 0, k +. 


This is the generalized biorthogonality condition of which we have already 
found special cases. 

Let us denote by D,; the first member of (18) and let us suppose that 
Dy, + 0 for each k. Then if we have a set of expansions of the form 


(19) fi(a, 8) = (z, 8), 1,2, «++, 
r=!) 


giving the simultaneous expansions of the n functions f,(2, s), 7 = 1, 2, 
-++, n, in terms of the y;“” (a, s), the coefficients c, may be found by mul- 
tiplying both sides of the equation 


aji(x, 8)2; (a, 8) + SEL t, s)2; (a, t) ldt, 
j=l 


summing as to 2 from 1 to n, integrating as to s from a to B and as to x from 
a to b, and simplifying by use of (18); thus we have 


n n 
Dryty = Sa" SEL 8) 
t=! j= 


+ fed t, s)2; (a, Ofi(a, s)dt|ds}da, v= 1, 2,3, 


t=1 


(20) 


It is clear that one may apply to (15) and (16) the Volterra limiting 
process and so replace these equations by others in which the discrete 
variable 7 has given way to an additional continuous variable, the resulting 
integro-differential equation containing a repeated integration, and that the 
formal expansion problem subsists in the new form. One might then employ 
a system of such equations and apply again to them the Volterra limiting 
process; and so on indefinitely. 

3. Differential and Integro-differential Problems Involving » Param- 
eters.—Let us consider the v systems of equations 


dyn; 
(1) dz j=1 
i= I, 2, 
a separate system being formed for each value k of the set 1, 2, 3, ---, v; 
and with these let us associate the v systems of adjoint equations 


(2) 2, (— Gongs — — — AVG, a= 1, 2, +++, me. 
j= 
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For fixed k and 2 multiply equation (1) member by member by z;; and equa- 
tion (2) member by member by y;:, add the two resulting equations member 
by member, sum as to 7 from 1 to n;, and integrate with respect to x from 
a to b; thus we have 


(3) + + + = 0, 


Let us now suppose that for each value of / we set up on the quantities 
yui(b), yxi(a), = 1, 2, mz, and separately on the quantities z;,,(b), 
Zri(a), = 1, 2, --+, nz, linear homogeneous boundary conditions (analogous 
to those developed in detail in § 1) of such sort that relations (3) are implied 
by the boundary conditions alone and that the y-problem and the z-problen 
have the same sets , Ao, -+-, A, of characteristic values, that the 
number of sets is enumerably infinite and that these distinct sets are 
denoted by the named symbols for s = 1, 2, 3, ---. Let the corresponding 
solutions be denoted by y;;“, 

In equation (1) let us replace the quantities y and \ by corresponding 
quantities y® and \®; in equation (2) let us replace the quantities z and 
by corresponding quantities 2“ and \“” where r + s. For a fixed k multiply 
the first resulting equation member by member by 2;;‘” and the second 
by yxi®, add the two resulting equations member by member, sum as to i 
from 1 to nz, integrate with respect to x from a to b, and simplify by what 
(3) becomes on replacing the quantities y and z by the corresponding quan- 
tities y and 2. In this way we are led to the following relations: 


i=! j=l 
Nk 


tal j=l 


a separate relation being formed for each value k of the set 1, 2, ---, ». 
Now the differences — ---, — A, are not simultaneously 
zero; hence the determinant of their coefficients in the system (4) must 
have the value zero, that is, we must have 


Nk 


So = 0, r+ 8, 


j=l 


where the term of the determinant explicitly written is that in the /th 
column and kth row. 

The last relation may be reduced to a more convenient form. In the 
kth row of the determinant let us replace the variable x of integration by 
t, and the variables 2 and j of summation by 2; and j; respectively. Then 
the relation reduces to 
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the determinant employed being that whose element in /th column and kth 
row is that explicitly written. Then if we employ the symbol D to represent 
the determinant 

(5) ivjy (E15 to, | (te) | 


whose element in /th column and kth row is that explicitly written, we have 
for the final form of our generalized biorthogonality conditions the following: 


ny Ny nv 
(6) t,=ljy=! 


x dt, = 0, r+ 8. 


In what follows we shall suppose that we have before us the case in which 
the first member of (6) is different from zero when r = s. 
Now let us consider the set of functions 


with regard to the question of their simultaneous expansion in the form 


the coefficients c, being the same for every function of the set. In case such 
an expansion exists we have a ready means of determining the coefficients. 
c,. The form of the first member of (6) suggests the operations to be 
performed. Let us denote by B,, the first member of (6) and let us sup- 
pose that B,, + 0 for r = 1, 2, 3, ---. Then for the determination of the 
c, we have the relations 
mm 
Bute = So 


A=! 
my Np 


(8) Diy (his to, t,) inte eee to, 7 t,) 


(tr) 22 (ta) 2°, (t, )} 
The expansion problem which is set by these formal results we may 
look upon as a generalization of that briefly treated at the beginning of § 2. 
From the fuller treatment of a related problem given in § 1 one sees readily 
how to proceed to a more precise formulation of the problem here briefly 
characterized. Moreover, it is clear that one may formulate the problem 
at once for a system of v equations of general orders and so include in one 
form (rather complicated, to be sure) the three differential equations 
problems already set. 
We saw in § 12 how the integro-differential equations problem (7), (8) 
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grew out of the differential equations problem (1), (2). In a similar way 
the problem (1), (2) of this section gives rise to a new problem associated 
with the equations 


Oyx(2, 8) = SF (2, 8, t) + A101, (2, 8, t) + 


(9) Ox 
+ d,a,x(a, 8, t)} yx(a, t)dt, 
(10) = SP {— aox(a, t, 8) — t, 8) 


— t, 8) t)dt 


With the guide afforded by the earlier part of this section and the special 

case of our present problem treated in § 2 one may proceed readily to set 

up the generalized biorthogonality conditions and by aid of them to deter- 

mine the coefficients in the formal expansion of a function of 2y variables. 
Corresponding to equation (6) we should thus have 


11 
A=l 


where the symbol D denotes the determinant of order vy whose element in 
ith column and fth row is If we denote the first member 
of this relation by C,, and suppose that C,, + 0 for every 7, then for the 


coefficients c, in the expansion 
(12) U2, °° $1, 82, ° = II (Xn, Sp) 
s=l h= 


we have formally 
(13) TT (Xn, th) dsidt, ds,dt,} 
+++ dz,, pom 1.3, 2, 


These results might, if it were desired, be also extended by means of the 
suggestion afforded by the problem associated with equations (15) and (16) 
of § 2 and the generalizations of it which are mentioned there. 

4. Problems Involving Difference and Integro-difference Equations.— 
Let us consider the adjoint systems of difference equations 


(1) Uu;(x + 1) Gy + + U;(2), a = 2, 
(2) 0;(2) = + + 1), a = 3, 2, 


where 6;; has the value 1 or 0 according as 7 is or is not equal to 7 and ¢ij 


i 


— 
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and y;; are functions of x which are analytic at infinity and vanish there 
to an order at least as high as the second. Under such hypotheses these 
equations have for a fixed value of \ fundamental systems* 


of solutions each function of which is analytic throughout the finite plane 
except for a set of singularities in a left half-plane; moreover, if x approaches 
infinity along any ray from the origin exclusive of the negative axis of reals 
or along any line proceeding to the right parallel to the axis of reals, we have 


lim u;;(7) = 6;;. lim = 6,;. 


-Any solution of either equation formed from the foregoing solutions of that 


equation by taking a linear combination of such solutions with constant 
coefficients will have the property that each function in it approaches a 
limiting value as x approaches infinity in the manner specified. We confine 
attention to such solutions of equations (1) and (2). 

Now if we multiply equation (1) member by member by v;(% + 1) and 
(2) by — u;(a), add the resulting equations member by member, and sum 


as to 2 from 1 to n, we have 
(3) DA{ui(x)vi(x)} = 0. 
i=l 


If the real part of a is sufficiently large we have wu; and v; analytic at every 
point « whose real part is not less than the real part of a. Hence in (3) 
we may sum as to x from a to infinity, where x runs over the values a, a + 1, 
a+ 2, --++; thus we have 


This is clearly analogous to equation (3) of § 2 and differs from it (apart 
from notation) only in having the limit © instead of the limit b. 

Let us now suppose that adjoint homogeneous linear boundary condi- 
tions are set up on the u;(), w;(a) on the one hand and on the v;(), v;(a) 
on the other hand so that relation (4) is implied by the boundary con- 
ditions and that this is done in such a way that the characteristic values 
of \ for (1) and the conditions on the w; are the same as those for (2) and 
the conditions on the v; and that the number of characteristic values is 
infinite. Let us denote the distinct characteristic values and corresponding 
solutions by A, ui, 15, k = 1, 2, 3, --> 

In (1) let us now replace \ and u; by A and u;, in (2) let us replace 
d and 2; by AM and v,;, where / + k; then we have two new systems. Let. 


* See AMERICAN JOURNAL OF MATHEMATICS, 35 (1913), 161-162. 
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us multiply the first resulting equation member by member by 2,” (x + 1) 
and the second by — u;“)(a), add the resulting equations member by mem- 
ber, sum as to 2 from 1 to n, sum as to x over the set a, a+ 1,a+ 2, --., 
simplify by aid of what (4) becomes when w; and »; are replaced by u; 
and »v, respectively, and divide by the non-zero quantity \ — \, 
Thus we have 


=0 j=! 


If we suppose that the quantity in the first member of (5) is (for all k) 
different from zero when / is replaced by k we shall have a ready means for 
the formal determination of the coefficients c;, in the (assumed) expansions 


(6) fix) = Neue), t= 
k=1 


the same coefficients ¢; being employed for each of the n functions. For 
this purpose we multiply both sides of equation (6) by y;;(x)o;(# + 1), 
sum as to 2 and 7 each from 1 to n, and sum as to 2 over the set a, a+ 1, 
a-+ 2, --+; employing (5) we come through readily to the relations 


n n 


= t=0 t=1 j=1 
(7) vila + + tj(a+1+t) 


Since the coefficients c;, are determined merely by the function-values 
at a discrete set of values of x it may seem at first sight that this is not 
likely to lead to an expansion theory of much interest. But it would be 
natural to confine attention to functions f;(7) having suitably restricted 
properties, especially in the neighborhood of infinity, and for such functions 
it may be anticipated that the theory is valuable. Moreover, it furnishes 
also a natural means of interpolation of a function defined initially only 
over the set of points a, a+ 1, a+ 2, ---, the function values at these 
points being the only ones needed in the determination of the coefficients 
c,. Again, it affords a natural means of formal expansion of functions 
defined only for positive integral values of the argument, such as is the 
case with the usual number-theoretic functions. 

Let us rapidly illustrate the foregoing theory by considering the problem 
associated with the equations 


(8) 1) — (1 ) ua) = 
©) (1-4 v(x) = 0, 


| 
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where pw? is the parameter. These equations have the solutions 


where I'(x) denotes the usual gamma function; and the general solution 
for each case is gotten by multiplying the particular solution by an arbitrary 
periodic function of period unity. But if u(x)[o(x)] is to approach a 
definite limiting value when x approaches infinity in the sense defined 
above this periodic multiplier must reduce to a constant; and such constant 
multiplier is clearly irrelevant to the problem in hand. Hence, for our 
purposes, we take the solutions given in (10). 

Now condition (4) reduces in the present case to u()v(%) — u(a)o(a) 
=(). As suitable boundary conditions implying this relation we may take 


(11) — = 0, 


where ¢ is a given non-zero constant. Now from (10) and the asymptotic 
properties of the gamma function we see that u(%o) = 1 = Then 
the condition of consistency for the u-problem and that for the v-problem 
both reduce to the following: 


(12) T(a — p)T(at p) = t{T(a)}?2. 


The roots of this equation give the common characteristic values of py? 
for the u-problem and the v-problem. [For each solution p of (12) there 
is also the solution — yu; but only their common square py? is a characteristic 
value for our problem. | 

For varying choice of a and ¢ the expansion problem which results 
brings a variety of properties to notice. One of the most interesting cases 
is that in which a= 3, —1=1/t=+1; in this case equation (12) 
readily becomes 


(13) cos T 


if one reduces by aid of the relation I'(s)I'\(1 — s) = z/sinas. Then the 
characteristic values yw? are readily determined in terms of a single one of 
them. Let yu be the smallest positive value of yu satisfying (13); then the 
characteristic values p? are 


w= 2n), n=0,41,42,43,---. 


Corresponding to these let us denote the solutions by u,(x), n(x) so that 
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we have* 


n = 0, 1, 2, 


The biorthogonality conditions now reduce to 


(15) ] un(k + 4) 


& = 0, Mm += Nn. 

Now we may apply to systems (1) and (2) a Volterra limiting process 
analogous to that by which we arrived at equations (7) and (8) in § 2, 
By so doing we are led to the integro-difference equations + 


(16) Au(a, s) = s, t) + AW(a, 8, t)} u(a, 
(17) — Av(z, = {o(a, t, 8) + AW(a, t, s)}o(a + 1, 
Corresponding to equations (4) to (7) we now have the following: 
(18) SE {u(%, s)v( s) — u(a, s)v(a, s)}ds = 0, 


(19) SSP fPy(atr, ts)m(a+r, k+l, 


r= 


(20) fle, = 8), 


+ r,t, s)\f(at+r, + 1+ 71, tdsdt 

+r, t, s)ux( s)v.(a+ 1+ 7, tdsdt 

k = 1, 2, 3, 

The foregoing work of this section is analogous to a part of that of § 2 
It is clear that the analogues of the other parts of § 2 exist here and can be 
readily developed. In fact, we can develop the theory of difference and 
integro-difference equations in forms analogous to all parts of the problems 
treated or suggested in §$1 to 3 for differential and integro-differential 
equations. 

5. Mixed Problems Involving One or More Parameters.—In view of the 
basic algebraic theory and the transcendental problems treated in §§ 1-4, 
it is clear that there must exist in the theory of integral equations expansion 
problems involving n parameters not only in the classic case when n = 1 
but also in the general case when n is any positive integer. It was my 
intention to exhibit the basic formule on which such an expansion theory 
may be based; but, while this paper was in course of preparation, an abstract 


* It should be observed that the negative values of n merely give us the same solutions 
over again. It is also worth noting that u,(x)/u:(z) is a rational function of x so that from 
the expansion of f(x) in terms of un(x) we have that of f(x)/ui(x) in terms of a certain 
interesting set of rational functions. 


t Here Af(z, s) = f(x + 1, s) — f(a, 8). 


eSS 
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of a paper by Mrs. A. J. Pell appeared in the Bulletin of the American 
Mathematical Society, Ser. 2, Vol. 26 (1920), p. 149, indicating the character 
of this expansion problem for the case of two parameters. Consequently 
we shall leave out the formule for this case. 

If we think of the several types of expansion problems—those for 
differential, difference, and integral equations—in intimate connection with 
the basic algebraic theory (developed in the memoir already referred to), 
it becomes apparent at once that the case of n parameters (for n > 1) is not 
confined to a set of n equations of the same sort. There is nothing to 
prevent one subset of the basic algebraic equations proceeding to differential 
equations as limiting forms, another to difference equations, another to 
integral equations, and so on. Thus we can see beforehand that we may 
formulate the expansion problem for a variety of mixed systems. In fact, 
we were led incidentally in § 2 to a special case of such mixed problems. 

Let us now indicate more precisely the nature of these mixed problems 


by considering the equations 
(1) dy: = >» (ai; + Abi; + + vdij)y;, 1 = 1, 2, 


dx 


m 


(2) As;(x) (ij + Api; + + V7 ij) 8; (2), i= 2; 
j= 
(3) u(x) = AS? E)u(E\dE + wf? L(x, + 
M (a, &)u(x, &)dé, 


associating with them their adjoint equations 


di 
(4) = (— aj; — — — 1= 1,2, 
j= 
™ 


(5) At;(x) = 24 (= Apji — vr 1), 1, 2, M, 


(6) = ASV K(E, + w L(E, + v M(E, x)v(E)dé. 
We suppose that the ¢, p, o, 7 functions are analytic at infinity and vanish 


there to at least the second order. 
From (1) and (4), from (2) and (5), by methods already frequently 


illustrated, we have 


(7) + yore + + Yntn = 0, 
(8) — si(a)ti(a)] = 0. 


We are now to set up homogeneous linear boundary conditions on y and on 
z so that (7) is implied by them; and on s and on ¢ so that (8) is implied 
by them. We shall suppose that this is done in such a way that the sets 
of characteristic values for the two problems [ (1), (2), (3) on the one hand 
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and (4), (5), (6) on the other ] are the same and that the number of the sets 
is enumerably infinite. Then we shall denote the corresponding values 


and solutions by 
AM, pH, 20, 4,0), yO, = 1, 2,3, 


In equations (1), (2), (3) let us replace the symbols for the parameters 
and the functions sought by the corresponding symbols with the super- 
script k; in equations (4), (5), (6) let us introduce similarly the distinct 
superscript /; and let us denote in order the six new equations so formed 
by (1), (2), -+°, (6). From (1), (4) and the boundary conditions implying 
(7) we have readily 


(9 = 


i=l j=l 


Similarly, from (2), (5) and the boundary conditions implying (8), we have 


(AM — XO) pita + w)s(at w)t(at1+w)+ - 


(10) 
+ (ph — yp) + w)si(a+ w)t(a+1+w)=0. 


w=0t=1 j=1 
From (3) and (6) we get likewise 


(11) AP — Sa KE, xu (x)v (E)dadé + - a 
+ vO) f? ME, (x)o (E)dadé = 0. 
Since the parametric differences \“ — ©, ---, involved in the last 
three equations are not simultaneously zero when k + / the determinant 
of the coefficients of these differences must have the value zero when k + 1. 
The relations thus arising express the generalized biorthogonality property 
of the solutions pertaining to the problem in hand. Once these are obtained 
the formal determination of the coefficients in the formal expansions 


(12) Fis(t1, = (a1)8; (as), 


of a set of mn functions of three variables is easily effected by the method 
already employed several times in this paper. 

The foregoing we have outlined briefly as an instance of a sort of thing 
which may be done in a great variety of ways. The formal aspects of the 
problem are easily developed whether we have a single equation of a given 
type, as in the problem here treated, or any number of equations of each 
type. Also, we may employ in the list of types integro-difference and 


n n n n 


ets 
ues 


“T'S 
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integro-differential equations. In view of the developments of the section 
following it will be seen that partial differential and partial difference equa- 
tions are also admissible. Similarly one may allow equations involving 
functions of two variables and the operation of differentiation with respect 
to one and that of differencing with respect to the other; and soon. In this 
vast totality most types are without special interest; but a few actually 
arise in problems which come naturally to our attention. 

Nor is this all. It is also possible to treat similarly equations in which 
functions are subjected to different operations with respect to the same 
variable. We illustrate this with one of the simplest cases. Let D; denote 
a set of n operations, the first k of them being identical with D and the 
remaining ones with A, where k is some one of the numbers 0, 1, ---, ”, 
and where D is the symbol for differentiation and A is the symbol for the 
operation Af(x) = f(a+ 1) — f(x). Then form the system of equations 


(13) Dyus(a) = (aij + 1, 2, +++, 0, 
j=l 


and associate with them the adjoint system 


(14) = (— — + 6), 1,2, ---, 0, 

j=l 
where e; is zero or unity according as D; is D or A. If we multiply (1) by 
v;(a + e;) and (2) by w;(x), add the resulting equations member by member, 
and sum as to 7 from 1 to n, we have a result which may be put in the form 


(15) = 0. 


Then, if appropriate conditions at infinity are satisfied, we have 


2, — u;(a)v;(a)} = 0. 
It is clear now how one may complete the formulation of the expansion 
problem; the work proceeds in the closest analogy with special cases of it 
in §§ 2 and 4. 

6. Problems Arising from Partial Differential Equations and Integro- 
differential Equations.—As typical of the general situation in regard to 
expansion problems arising from the theory of partial differential equations, 
let us consider that which is to be associated with the equation 


(1) L(u) + AL;(u) = 0 
and its adjoint 
M(v) + AM,(r) = 0, 


n 
n 
n 
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where 
Oxdoy 


lie) ma + 2b 


(i) 


M(v) = (ar) +2 


M,(v) = + ru, 


where the symbols wu, v, a, b, c, a, B, Y, p, 9, r denote functions of the in- 
dependent variables x and y. 

Let us multiply both members of equation (1) by » and of equation (2) 
by — wand add the resulting equations member by member; thus we have 


0 Ou Ov Ou Ov da Ob 
=| use) +0( )+ (a+ | 


du av Ou Ov 


Let S be any closed region in the xy-plane whose boundary is cut twice 
and only twice by a line through it and parallel to either axis. Then if 
we integrate over S both members of equation (3) we have a relation 
which reduces readily to the form 


Ou Ov. Ou Ov Oa ab 
Ou Ou Ov 


Y=Cz 
+ — | dx = 0, 
Oy Ox 


where r, and s, are limiting values of x in S for a given y and r and s are 
the extreme values of x and where p;, oz, p, o have like meanings under 
interchange of réles of 2 and y. 

If in the first and second integrals in (4) we introduce the variable ¢ 
in place of the given variable of integration by means of the relations 


respectively, each integral reduces to an integral with respect to ¢ from 


dc Ob 
+({B+)q By | 0. 
s—r 
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0 to 1, and hence (4) itself may be put in the form 
(5) JS o'Bdt = 


where B is a sort of bilinear form in the functions u and v and their first 
derivatives, the arguments of the functions being specialized in the way 
implied by the transformations involved (so that each of them is a function 
of a single variable £). 

We now require linear homogeneous boundary conditions on u and v 
separately, independent of A, such that equation (5) is implied by the 
boundary conditions alone. We shall presently make clear the nature of 
these conditions by the explicit treatment of a particular case. For the 
present we shall suppose that these conditions have been set up in such 
way that the w-problem and the v-problem have the same characteristic 
values \ and that these are enumerably infinite in number. This done, we 
denote the characteristic values \ and the corresponding solutions by \,, 
U;, = 1, 2, 3, 

From (1) and (2) we now have the relations 


+ (u;) = 0, 
M(v;) + (0;) = Ai — Aj) 


If we multiply the first of these equations member by member by 2; and 
the second by — w;, add the resulting equations member by member, and 
then integrate over S, the resulting first member has the value zero on 
account of the boundary conditions, and hence the second member also 
has the value zero. Omitting the non-zero factor \; — \;, this gives us 
the relation 


(6) S (0;)dady = 0, 0 

Similarly, we have 

(7) S S = 0, 1 #7. 


These two (equivalent) sets of relations express the generalized biorthogo- 
nality conditions for the problem in consideration. 

If we suppose that the integral {” f-u:M,(v;)drdy taken over S is dif- 
ferent from zero for every i we have a ready means of determining formally 


the coefficients c;, in the formal expansion 


8) fle, y) = 


of a function f(z, y) of two variables in terms of the functions u;. In fact, 


in- 

2) 

ve 

0. 

if 

n 

> 
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we obtain readily the relations 


SSG 


SS wile, My(o,)dedy 


In order both to make the method of procedure more explicit and to 
illustrate the greater simplicity of the self-adjoint case of the foregoing 
problem let us consider an instance associated with a class of equations 
of great importance in theoretical physics, namely, equations of the form 


(10) A(g + hju = 0, 


where g is a function of 2 alone and h is a function of y alone. For the 
region S we take now the square 0 = 2 =1,0 Sy =1; and we use the 
boundary conditions 


(11) u(t, 1) — u(t, 0) = 0, u(1, t) — u(0, t) = 0, 
u,(t, 1) — u,(t, 0) = 0, u,(1, t) — u,(0, t) = 0, 


where the subscript y or x denotes the partial derivative with respect to 
this variable. It is easy to verify that problem (10), (11) is self-adjoint; 
that is to say, if one writes the same problem in v one finds that these two 
problems, one in v and one in wu, satisfy those conditions already imposed 
in the more general problems associated with (1) and (2). 

Equation (6) now has the special form 


(12) Sc S dg + hjuwdrdy = 0, 1 +49, 
so that for the coefficients c; in (8) we have the simpler values 


SV + h){ux(x, y) }Pdady 


Instead of considering all the solutions of (10) which satisfy conditions 
(11) we may fix our attention on any convenient class of them, as for in- 
stance those solutions which can be written in the form 


(13) 


(14) u(x, y) = p(x)o(y), 


a class of particular solutions of great importance in theoretical and applied 
mechanics. Substituting u(x, y) of the form (14) into (10) and (11) we 
find that a constant p exists such that 


(15) p’’(x) + Ag(x) p(x) + wp(x) = 0, p(1) — p(0) = 9, p’(1) — p’(0) = 9, 
(16) o’(y) + Ah(y)o(y) — wo(y) = 0, o(1) — (0) = 0, o’(1) — o’(0) = 0. 


he 
he 
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Each of the last two differential systems is self-adjoint. Together they 
form [if written in a single variable 2] a system of two differential equa- 
tions (with boundary conditions) in two parameters \ and yw. General 
problems of this type have already been treated in $3. It was through 
the emergence of the multiple parameter problem in connection with partial 
differential equations that I was first led to seek the formulation of the 
general multiple parameter problem for algebraic systems and also for 
the transcendental cases treated in this memoir. 

Let us derive the equations for the characteristic values of \ and yp in 
equations (15) and (16). Let pi(v, A, w) and pe(a, A, w) be a fundamental 
system of solutions of the differential equation in (15) and seek a solution 
p(x, A, wu) so that the boundary conditions shall be satisfied. Constants 
c, and ¢, are to exist such that 


p(x, = C1pi (2, A, + Copo(2, 
hence we have through use of the boundary conditions in (15) the relations 


C1{p1 (1, (0, + C2 { (i, r, =" (0, = 0, 


The condition of consistency of these equations is that 


(1, A, px (0, p2'(1, A, p2'(0, A, 
a condition on and yw alone. Similarly, if oi(y, uw) and w) are 
a fundamental system of solutions of the differential equation in (16), 
we are led to the following second condition on X and u: 


o1/(1, d, w) — A, A, — A, 


(18) 


The common solutions of equations (17) and (18) are the characteristic 
values for problem (15), (16). They are also the characteristic values for 
problem (10), (11) when the solutions w are restricted to be of the form given 
in (14). 

Let us consider for a moment the special case in which g = h = 1. 
Then equations (15) and (16) have solutions (among others) of the form 


p = e*, o = 


where @ and 8 are constants. The boundary conditions then require that 
we must have 


e—1=0, a(e*— 1) = 0, —1=0, — 1) = 0. 


Hence a = 2mmi and 8 = 2nzi where m and n are integers, positive or 


to 
hig 
| 
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negative or zero. Then corresponding characteristic values \ and p [per- 


haps not the complete set | are given by the relations \ + u = — 4m?z? and 
A— w= — 4n’r’, or 
A= — + n’), — — n’), 


where m and n are integers. The corresponding solutions u(x, y) may now 


be written 


The expansion of f(a, y) may be written 
+x 


f(z, y) = iy, 


m,n=—® 


For a real-valued function f(x, y) it is easy to transform this into an expan- 
sion in sines and cosines of integral multiples of x and y with real coefficients. 

Let us consider the extension of our problem to adjoint systems of partial 
differential equations involving a single parameter. As a typical instance 
we take the system 


Ou; , < Ou; du; 
(52+ bi; + cus) 


Ou; Ou; 


j=l 


and its adjoint system 


Let us multiply the first of these by v; and the second by — u;, add the 
resulting equations member by member, and sum as to? from 1 to n; thus 


we have 


2.10 Ou; Ov; Ou; Ov; 


oo 
+ + — + Uri | = O. 

After the method employed in the first part of this section we may now 
obtain from (21) the condition which must be realized as a consequence of 
linear homogeneous boundary conditions on the u; and the 2; separately. 
In (21) we integrate over a two-dimensional region S of the xy-plane, per- 
forming the integration with respect to x for those terms affected by the 


— 

077; , 079; < Ov; 00; Oa;;  0b;; 

(20) ay ay +(¢ dx 


nd 


OW 


we 
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operator 0/dx and with respect to y for those affected by the operator 0/dy. 
In this way we have a relation analogous to (4) above. Just as in the 
preceding case we introduce a variable ¢ so as to reduce the two integrals 
to a single integral of the form 


(22) Jo Bdt = 0 


where B is a sort of bilinear form. The boundary conditions must be set 
up so as to realize this relation; they may be set up so as to realize the more 
special relation B = 0. 

Without attempting to characterize these boundary conditions in general 
[suggestions for forming them are afforded by the earlier part of this 
section and by the treatment of a like problem in § 1] we suppose that 
they have been chosen so that those on the wu; are linear and homogeneous 
in the uw; and their first derivatives and so that those on the 2; are linear and 
homogeneous in the 2; and their first derivatives and that this has been 
done in such way as to realize condition (22) as a consequence of the bound- 
ary conditions. We shall suppose also that the characteristic values \j, 
Az, Az, Of the w-problem are the same as those for the v-problem and 
that these values are infinite in number. The solutions corresponding to 
we denote by u;“, o;. 

In (19) let us write A, and wu; for \ and w;; in (20) let us write \; and 
vi for X and v;, where k +1. Multiplying the first of these resulting 
equations by »;‘” and the second by — u;“, adding member by member, 
summing as to 2 from 1 to n, integrating as to x and y over the region S, 
simplifying by use of the relation to which (22) reduces when w; is replaced 
by wu; and 2; by v;“ [a relation which is valid in view of boundary conditions 
of the sort described ], we have a relation which becomes the following 
when the non-zero factor A, — dz is dropped. 


i=l j=l Ox 


Ox Oy 


(23) 


Similarly we have the equivalent relation 


(24) SS + + | dedy = 0, k + 1. 


1 j=l 


Now let fi(z, y), 1 = 1, 2, ---, n, be a system of n functions of the two 
variables x and y and let us seek simultaneous expansions of the form 


v=1 
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where the c, are independent of 7 as well as of x and y. We seek to deter- 
mine the c, on the supposition that the expression in the first member of 
(24) is different from zero whenever k = 1. We multiply both members 
of (25) by 
Ox Oy 

sum both members as to 7 from 1 to n a" as to 72 from 1 to n, and integrate 
in the resulting equation as to x and y over the region S._ In view of (24) 
it is seen that c¢; is the only coefficient left in the second member and that 
we have for c; the value 


n n x ra] J 


i=l j=l 


v, y)dady 


k= 1,2,- 


In case a;; = 0 = £6,; relations (24) and (26) reduce to the notably 
simpler forms 


Ck = k= I, 2, 3, 


i=l j= 


Just as the problems in ordinary differential equations with one param- 
eter can be generalized, as we have seen, to the case of any finite number of 
parameters so the problems treated in the foregoing part of this section 
can be extended to partial differential equations and systems in any finite 
number of parameters. Moreover, any finite number of independent 
variables may be employed in place of the two variables x and y of the fore- 
going treatment. 

If we use the Volterra limiting process for passing from a finite number 
of differential equations to an integro-differential equation we are led from 
equations (19) and (20) above to certain adjoint partial integro-differential 
equations. For the sake of simplicity we shall confine our attention to the 
case in which the functions a;;, b;;, a;:;, 8;; are all identically zero. Then the 
integro-differential equations are the following: 


(27) 8) _ SP y, 8, t) + y, 8, t) Ju(a, y, tdt, 
2 


2 
(28) +— = fl e(a, y, + y, t, 8) Jo(a, y, t) 


of 
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These equations may be treated by processes closely analogous to those 
by which (19) and (20) were treated. We summarize the main results. 
Corresponding to (23) and (24) we now have 


(29) SP y, t, (a, y, thu (a, y, s)dsdtldady=0, k=l. 
Then if we seek expansions of the form 


(30) f(x, y, 8) = 2, cw (x, y, 8) 
v= 
we have for the coefficients c; the values 


SITIE SEV, y, 8, thu”(x, y, Of(a, y, s)dsdt |dady 


7. On a Problem Arising in the Theory of Vibrating Plates.—An impor- 
tant differential equation in the theory of vibrating plates is the following: 

Ax’ dy? 


Here a is a constant and é is a function of x, y, t. In the theory of vibrating 
plates one desires to have a solution £(2, y, t) of (1) which for ¢ = 0 reduces: 
to a given function f(a, y) so that we have &(a, y, 0) = f(z, y). Such solu- 
tions &(a, y, t) have been found in a variety of ways in terms of particular 
solutions of the form 


where n is a constant. If we substitute this value of £ into (1) we have for 


u the equation 


=, dtu 
3) AU = Au = 0, 
where \ = — n?a. We shall formulate the expansion problem as con- 


nected with equation (3) in the general form* suggested by the standard 
procedure which we follow under the guidance of the basic algebraic theory 
which directs our investigation. 

The equation adjoint to the equation L(u) + Aw = 0 is obviously the 
equation L(v) + A(v) = 0. We have the following readily verified basic 
identity : 

* A number of special cases of this problem have been treated (with important results) 
in connection with the usual theory of vibrating plates. 
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vL(u) — uL(v) = 


E _ 9° _ dvdu , 
2 3 
2 3 


— su 
In this identity s denotes an arbitrary parameter. 

Let S be any closed region of the zy-plane whose boundary is cut twice 
and only twice* by a line parallel to either axis and passing through the 
region. Then if u and v are two solutions of (3) we have 

| 
0= SSf[oL(u) — uL(v)] dxdy = Sf — Bit Bz ; drdy, 
S S dx dy J 


where B,; and By, denote the first and second brackets expressions in the 
second member of (4). From this we have readily the relation 


(5) (Bi + Bs = 0 
where r,, Sy, 7, 8, Pz, Tz, p, © have the same meanings as in connection with 


equation (4) of §6. By the same transformation as is employed in § 6 
we can reduce relation (5) to the form 


(6) Si Bdt = 0 


where B is a sort in bilinear form in the functions wu and v and their deriva- 
tives up to the third, the arguments of the functions being specialized in the 
way implied by the transformations involved (so that each of the functions 
is a function of the single variable f). 

We now require (as usual) linear homogeneous boundary conditions on 
u and v separately such that equation (6) is implied by the boundary condi- 
tions alone. A particular class of cases of great importance is that in which 
the boundary conditions imply the stronger relation B= 0. Leaving off 
for the moment the actual formulation of these boundary conditions let 
us suppose that they have been set up in such a way that the problem for 
the equation L(u) + Au = 0 and the w-conditions has the same character- 
istic values \ as the problem for the equation L(v) + Av = 0 and the 2- 
conditions and that these characteristic values are enumerably infinite in 


* The reader will see what modifications of the treatment are sufficient to deal with 
the case when these lines cut the boundary any finite number of times. 


ice 


he 


CARMICHAEL: Boundary Value and Expansion Problems. 267 


number. We denote the characteristic values \ and the corresponding 

solutions by Ai, wi, ;, 7 = 1, 2,3, -++. Since equation (3) is self-adjoint it 

is seen that wu; = v; for every 7 provided that the w-conditions and the v- 

conditions are the same. In this case the entire u-problem is self-adjoint. 
Now we have the relations 


L(u;) + = 0, L(v;) + A; = Qi — 


Let us multiply the first of these equations by v; and the second by — u;, 
add the resulting equations member by member, and then integrate over S; 
the resulting first member has the value zero on account of the boundary 
conditions, and hence the second member also has the value zero. Omitting 
the non-zero factor A; — dj, this gives us the relation 


(7) SS upjdedy = 0, #). 


S 


For the self-adjoint case this reduces to 
(71) SS uujdedy = 0, +). 


If we suppose the integral in the first member of (7) to be different from 
zero when j = 7 (as it certainly is in the self-adjoint case) we have a ready 
means of determining the coefficients c;, in the formal expansion 


(8) f(a, y) = y) 
of a function f(x, y) of two variables in terms of the functions uz. In fact 


we obtain readily the relations 


SS f(x, y)ou(x, y)dady 
(9) SS una, ydady 
S 


In the self-adjoint case we have the simpler expressions 


SS, f(x, y)dady 


k = 1,2, 


Suppose that we desire a solution £(2, y, ¢) of (1) which for ¢ = 0 reduces 
to f(a, y) so that we have formally 


E(a, y, 0) = cxun(a, y), 
k=l 


where the c;, have the values just obtained. Such a solution is formally 


i! 
t 
| 
| 
1e 
h 
| 
| 
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the following: 
E(a, y, t) = cxun(a, yev 
This is the actual solution desired provided that the boundary conditions 
are such as to render the solution unique and provided that the requisite 
conditions of convergence are satisfied. We have thus completed the 
formulation of the expansion problem for the present case. 

It is desirable to observe how easily one can extend the formal theory 
associated with equations (3) to (9) inclusive to that depending on the 
equation L(u) + Agu = 0 in place of equation (3), where g is a function of 
x and y. The boundary conditions will be set up exactly as before, the 
function g not intervening in this part of the discussion. In the paragraph 
containing (7) we shall have a slight modification of formulz and there 
will result the generalized biorthogonality conditions 


SS = 0, + 3; 


S 


and these give rise to the relations 
SS gla, yf (a, ydady 
S 

SS g(x, yur(a, 
S 


Ck 
for determining the coefficients ¢; in the formal expansion (8). 
It is clear that the formal expansion theory connected with the equation 
L(u) + Agu = 0 can be readily carried over to the system 


n 
= (au + rgiu), i= 1,2, +++, 0, 
j=l 


and to the simultaneous expansion of n functions in terms of the solutions 
obtained. This in turn can be made to yield heuristically the formulation 
of a corresponding problem in partial integro-differential equations. These 
two problems in one parameter can thence be easily extended to problems 
in any finite number of parameters. The method of making these exten- 
sions is sufficiently obvious in view of the methods already employed in this 
paper in closely related problems. 

In order to make clear the nature of the boundary conditions associated 
with equation (3) and described above through the intervention of equation 
(6), let us set them up explicitly and generally for the case when S is a 
rectangle and relation (6) is to be satisfied by requiring that B shall be 
identically equal to zero. We define S as that rectangle in which a = x = }, 


<= 


| 
| 
t 
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Let us now write more explicitly 
(10) Bile, y), Be(x, 9), 


respectively, for the forms in the first and second brackets expressions in 
relation (4). Then for B we now have the value 


B = (d — c){Bi(b, e+ (d — c)t) — Bila, e+ (d— e)t)} 
(11) + (b— a){B.(a+ (b — a)t, d) — Bo(a + (b — abt, c)}. 


Now let us think of B,(2, y) as a sort of bilinear form in the two sets of 
variables 


and 2,- 


Ox Oy? Ox? Ox dy? Oxdy" 


Let us form the matrix | 
of the ith column and jth row is the coefficient of the product of the ith 
function in the first set by the jth function in the second set. Then we have 


‘a;;|| of this bilinear form where the element a;; 


0 0 0 s8 

0 0 -1 —s 0 0 

0 0 0 0 
sl 0 2-s 0 00 0 
0 0 0 0 0 

s-2 0 0 00 0 


It is easy to show that this matrix is of rank 4; for the last row may be 
made to consist of elements zero by adding to it (s — 2) times the second 
last row, the fourth row may be made to consist of elements zero by adding 
to it (s — 2) times the third row, the last column may be made to consist 
of elements zero by subtracting from it s times the second last column, and 
the fourth column may be made to consist of elements zero by subtracting 
from it s times the third column, the matrix so obtained being the following 
(which obviously is of rank 4): 

0 0 0 0 1 

00-1 0 0 OF 

01 000 

0 0 000 


—-1 0 0000 
00 000 0) 


In By(a, y) we have the six variables wv and the six variables v listed in 
(12). In B.(a, y) we have similar six variables wu and six variables v. The 
matrix of the bilinear form B2(x, y) is of rank 4, as one may readily verify 
by the method used for B,(x, y). Hence we may look upon B, defined as 
in (11), as a bilinear function of 24 variables uw and 24 variables v, and this 


| 
e 
| 
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bilinear form is of rank 16. Hence B may be put into the form 


16 


(13) B= U;(u) V(r), 


where for each 2 U;(w) is a linear homogeneous function of the 24 variables 
u and V;(v) is a linear homogeneous function of the 24 variables v. The 
set U;(u), i = 1, 2, ---, 16, is linearly independent, and so also is the set 
Vi(v), 1= 1, 2, ---, 16. Then the boundary conditions which we take 
for the u-problem are 


(14) U,(u) = 0, 1,2, 8, 
while those for the v-problem are 
(15) Viv) = 0, «= 9, 10, ---, 16. 


[For some parts of the theory one may take any number & of conditions 
U,(u) = 0 not greater than 15, 2 = 1, 2, ---, k, and then use for the v- 
conditions V;(v) = 0,2 = k+1,---,16. But the form of conditions given 
in (14) and (15) is needful for certain parts of the development. | If the 
conditions on v in (15) are the same as those on wu in (14) the problem is 
self-adjoint and the maximum of elegance is secured. 


UNIVERSITY OF ILLINOIS, 
October, 1920. 
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RECIPROCITY IN A PROBLEM OF RELATIVE MAXIMA AND 
MINIMA.* 


By James K. WHITTEMORE. 


Introduction.—Most students of the differential calculus have doubtless 
observed that the two following problems have the same solution: To 
determine the shape of a rectangle of given perimeter and maximum area; 
To determine the shape of a rectangle of given area and minimum perimeter. 
While in this and similar cases a quite elementary explanation may be given 
of the exchange of maximum and minimum corresponding to the exchange 


" of area and perimeter it seems of interest to consider the analytical problem 


of maximum or minimum value of a function of two variables, subject to 
the condition that a second function of the two variables have a constant 
value, and to determine when the exchange of the réles of the two functions 
results in the exchange of maximum and minimum. In this paper this 
question of reciprocity is discussed for the simplest problem, two functions 
of two variables. It is suggested that a similar discussion of the more 
general problem with more than two variables and with two or more than 
two functions might give results of considerable interest. The methods 
are analytical but are based largely on geometrical intuition. It is intended 
that the discussion should be complete, that is that all necessary and suffi- 
cient conditions should be given, in so far as this is possible with the use of 
the first and second partial derivatives of the two functions. In the first 
section the general case is considered, simple necessary and _ sufficient 
conditions for a relative extreme, maximum or minimum, obtained and 
stated in both analytical and geometrical form, and the reciprocity condition 
determined. In the second and third sections the more complicated excep- 
tional cases are discussed. The fourth section is devoted to a study of the 
invariant properties of the conditions previously established. In the last 
section we give examples illustrating the theory developed in the preceding 
sections. 


/ 


$1. The General Case.—Consider two functions ¢ and y of the two 
real variables, x and y; suppose that in the neighborhood of P(20, yo) 
both ¢ and W are single valued, real, and have continuous partial derivatives 
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bilinear form is of rank 16. Hence B may be put into the form 
16 

(13) B= U,(u) V(r), 


where for each 7 U;(u) is a linear homogeneous function of the 24 variables 
u and V;,(v) is a linear homogeneous function of the 24 variables v. The 
set U;(u), 7 = 1, 2, ---, 16, is linearly independent, and so also is the set 


Vi(v), = 1, 2, ---, 16. Then the boundary conditions which we take 


for the u-problem are 


(14) Uu) = 0, 1=1, 2, 8, 

while those for the v-problem are 

(15) Viv) = 0, = 9, 10, ---, 16. 

[ For some parts of the theory one may take any number k of conditions 
U;(u) = 0 not greater than 15, 2 = 1, 2, ---, &, and then use for the v- 


conditions V;(v) = 0,4 = k+ 1, ---,16. But the form of conditions given 
in (14) and (15) is needful for certain parts of the development. | If the 
conditions on v in (15) are the same as those on wu in (14) the problem is 
self-adjoint and the maximum of elegance is secured. 


UNIVERSITY OF ILLINOIS, 
October, 1920. 
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RECIPROCITY IN A PROBLEM OF RELATIVE MAXIMA AND 
MINIMA.* 


By James K. WHITTEMORE. 


Introduction.—Most students of the differential calculus have doubtless 
observed that the two following problems have the same solution: To 
determine the shape of a rectangle of given perimeter and maximum area; 
To determine the shape of a rectangle of given area and minimum perimeter. 
While in this and similar cases a quite elementary explanation may be given 
of the exchange of maximum and minimum corresponding to the exchange 


" of area and perimeter it seems of interest to consider the analytical problem 


of maximum or minimum value of a function of two variables, subject to 
the condition that a second function of the two variables have a constant 
value, and to determine when the exchange of the réles of the two functions 
results in the exchange of maximum and minimum. In this paper this 
question of reciprocity is discussed for the simplest problem, two functions 
of two variables. It is suggested that a similar discussion of the more 
general problem with more than two variables and with two or more than 
two functions might give results of considerable interest. The methods 
are analytical but are based largely on geometrical intuition. It is intended 
that the discussion should be complete, that is that all necessary and suffi- 
cient conditions should be given, in so far as this is possible with the use of 
the first and second partial derivatives of the two functions. In the first 
section the general case is considered, simple necessary and sufficient 
conditions for a relative extreme, maximum or minimum, obtained and 
stated in both analytical and geometrical form, and the reciprocity condition 
determined. In the second and third sections the more complicated excep- 
tional cases are discussed. The fourth section is devoted to a study of the 
invariant properties of the conditions previously established. In the last 
section we give examples illustrating the theory developed in the preceding 
sections. 

§ 1. The General Case.—Consider two functions ¢ and y of the two 
real variables, x and y; suppose that in the neighborhood of P(2, yo) 
both ¢ and y are single valued, real, and have continuous partial derivatives 
of the first and second orders. We write 


=»% 
* Presented to the American Mathematical Society at New York, April 24, 1920. 
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and denote the values of these functions at P by wo and %. We consider 
the two problems: 

A. To determine when wp is an extreme of ¢ subject to the condition, 
y = %. 

B. To determine when % is an extreme of w subject to the condition, 
go = Uo. 

It is necessary to specify exactly what is meant by the statement, 
Uo is, for example, a maximum of ¢ subject to the condition ¥ = v. This 
statement shall in the following pages be understood to mean: (1) There 
exists a continuous set of real values x, y, of which set 20, yo is an interior 
point, and of which all points satisfy the equation Y = v; (2) For every 
point of this set, different from 20, yo, uo > g(x, y). Consider with this 
definition of relative maximum the example 


e=y, warty, 
Xo = Yo = U = % = OV. 


The condition, 2? + y? = 0, gives g = — 2’, and 0 is a maximum value of 
y considered as a function. of the independent real variable 2, but according 
to the definition above is not a maximum relative to the condition, 
x’? + y? = 0, for there is not a continuous set of real values, x, y, including 
0, 0, satisfying the equation of condition. In some circumstances, as will 
appear in § 2, there are two or more sets of real values 2, y satisfying the 
condition.* It may happen that wo is a maximum on both sets, supposing 
that there are two such sets, a maximum on one and not on the second, or 
not a maximum on either. 

Evidently problems A and B are interchanged by exchanging ¢ and y. 
Problem A may be stated in geometrical form as follows: We regard x and 
y as ordinary rectangular coérdinates in a plane which we find it con- 
venient to call horizontal; the problem is that of an extreme of the ordinate 
u of the curve of intersection of the surface, S, g = u, with the cylinder, 
y =. Suppose first that P is not a singular point of the curve C in the 
ay plane, Y = v; then one branch of C passes through P, and there is no 
essential restriction in supposing that, for P, Y, + 0. Evidently a neces- 
sary condition for an extreme in A is that the tangent to the curve of inter- 
section of S and cylinder at Q (ao, yo, Uo) be horizontal. For this curve 


* Usually in the following pages we consider only such a set of points (x, y) as form a 
curve through P having at P a continuously turning tangent. 
t See § 4. 
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It is necessary that, at P, 
J = = 0. 


The “critical points,” such as 2, yo, are determined in the usual elementary 
solution of A by solving simultaneously y = v) and J = 0. 
For P we have 


du 
Ly + Puy 2y Pry Pray y” = 
dz dx 

dx? Wy + + “y Wry + Wer = 


If, at P, g, = 0 it follows from J = 0 that ¢, = 0. This case we consider 
in § 2, supposing now ¢, + 0. We have at P 


Wy i Wy \ Py 


— + Cre = Vyy — Web + | 
Gy vy 


since J = 0. We may write for P 


Pu 


where y,” and y,” are the values at P of d’y/dx* for the curves in the zy 
plane, ¢ = uw and Y = w respectively. For problem B we should evidently 
have at P the necessary condition J = 0, and 

— 44”) 

We may state the following: With the hypothesis that each of the real 
curves in the zy plane, g = wp and w = %, has an ordinary point at P, a 
condition necessary in both A and B is that these curves be tangent at P, 
J = 0; if this is the case a sufficient condition in both A and B is that the 
curves do not osculate at P; if both the necessary and sufficient conditions 
are satisfied A and B have like or unlike extremes as ¢, and y, have different 
signs or the same sign at P. If the two curves osculate at P, so that the 
sufficient condition fails, problems A and B cannot be completely discussed 
by use of the first and second derivatives; we shall say in such a case the 
discussion fails.* 

§ 2. Exceptional Cases.—We consider in this section the cases of prob- 
lem A excluded in the preceding discussion. The conditions stated in each 


* For example see § 5. 
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case are for the values 20, yo. 
I. Ye = yy = 0, ¢ and ¢, not both zero. 
Il. g. = gy = 0, ¥z and y, not both zero. 

Ill. = gy = = fy = 0. 

To put the problem in geometrical form we consider 2, y as before to 
be rectangular coérdinates in a horizontal plane, and consider the inter- 
section of the surface S with the cylinder, y = v, a cylinder with vertical 
elements erected on the curve, C, Y = v, in the xy plane. This curve has in 
I and III a singular point at P. We assume that not all three of the second 
partial derivatives of y vanish at P, so that this point is a double point. 
This point is an isolated point with imaginary tangents, a cusp or osculating 
point or isolated point with one real tangent, or a point of intersection of two 
real branches of the curve with distinct real tangents, as AY = W2,? — Wr2Wry 
is negative, zero, or positive at P. The slopes of the tangents are the values 
of y’ satisfying 

+ + Ver = 0. 


In II and III the tangent plane to S is horizontal at Q. The total curvature 
of the surface at this point is positive, zero, or negative, for all cases, as Ag 
is negative, zero, or positive, and the slopes of the horizontal projections of 
the asymptotic tangents at the point are in all cases given by the values 
of y’ satisfying 

Sy + 2y'Gry + = 0. 


In accordance with the limitations set for our discussion we suppose that 
not all of the second partial derivatives of ¢ vanish at Q, or, in geometrical 
terms, that @ is not a flat point of the surface S. 

Case IL—y, = Py = 0, g, and ¢g, not both zero. P is a double point of 
C; the tangent plane to the surface S at Q is not horizontal, but contains a 
single horizontal direction given by y’¢, + ¢: = 0. We assume, as in the 
first section, g, + 0. 

Ta. AY < 0. No real branch of C passes through P, and wu is not an 
extreme of ¢. 

Ib. AY = 0. Suppose the point P is a cusp or an osculating point of C. 
If P is a cusp and if the cuspidal tangent is not parallel to the horizontal 
tangent to S, at Q, that is if 


D ge Vyy — + Wax + 0, 


the ordinate wu of a point moving on the curve of intersection of surface and 
cylinder will increase to uw then decrease, or decrease to wp then increase, 
since du/dx + 0 at P, and up is an extreme of ¢. The discussion fails in this 
case to distinguish between maximum and minimum. [If P is an osculating 
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point of C and if a path with no cusp at P be followed a necessary condi- 
tion for an extreme is Dyy = 0. The discussion fails for further conditions, 
and fails also to distinguish between cusp, osculating point, and isolated 
point. 

Ice. AY > 0. The point P is a double point of C with two real distinct 
tangents. In this case we consider only a path with a continuously turning 
tangent, that is a path projected into a branch of C. Evidently, for one 
branch of C, uo will not be an extreme of ¢. In order that wp shall be an 
extreme for one branch of C it is necessary that Dyg = 0. Further dis- 
cussion of this case fails. 

Case II.—g, = gy, = 0, wz and y, not both zero. The point P is a 
simple point of C; the tangent plane to S at Q is horizontal. In all cases 
under II, since, at P, y,” is finite with the assumption y, + 0, 


Pu _ Ds, 
de 


Ila. Ag < 0. The total curvature of S is positive at Q, and uw is an 
absolute and consequently a relative extreme of y, a maximum or minimum 
as D,w + 0 is negative or positive. It is evident that ¢,, and gy have 
the same sign as D,y. 

IIb. Ag = 0. The total curvature of S vanishes at Q, which is a para- 
bolic point, since we have excluded the case of a flat point. There is at Q 
a single asymptotic tangent given by 


+ 2y' Gry + = 0. 


It is clear geometrically that wo is a relative extreme of ¢ if the tangent to 
C at P is not parallel to the asymptotic tangent, that is, if Dsy + 0. Asin 
IIa the nature of the extreme is given by the sign of D,y, the same as that 
of ¢zz and gy, when neither of the latter vanishes. If D,y = 0 the discussion 
fails. 

IIc. Ag > 0. The total curvature of S is negative at Q, which is accord- 
ingly a hyperbolic point. There are at Q two distinct real asymptotic 
tangents given by the same equation as in IIb. If the tangent to C at P 
is not parallel to either asymptotic tangent wp is an extreme of ¢, whose 
nature is given as in Ila and b by the sign of D,y, but not by the sign of 
Grr OF Gy. If Dyy = 0 the tangent at P is parallel to one of the asymptotic 
tangents and further discussion fails. 

Case III.—o, = gy = vz = Wy = 0. 

The point P is a double point of C; the tangent plane to S at Q is hori- 
zontal. 

IIIa. Ag < 0. The total curvature of S is positive at Q. 
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1. AY <.0.. No real branch of C passes through P, and wo is not an 


extreme of ¢. 


2. AY = 0, 3. AY > 0. Since w is an absolute extreme of ¢ it is a 
relative extreme for any path, a maximum or minimum as ¢,, is negative or 
positive, if P is not an isolated point of C in 2. 

IIIb. Ag = 0. The total curvature of S vanishes at Q, and there is a 
single asymptotic tangent. 

1. AY <0. There is no extreme of ¢ at w. 

2. AY = 0. The curve C has a single tangent at the double point P. 
It is evident geometrically that wo is an extreme of ¢ if P is not an isolated 
point of C and if the tangent at P is not parallel to the asymptotic tangent 
at Q. 

To express the last condition analytically consider the values at P(2o, yo) 
of the two polynomials in y’, 


yo + 2y’ Gry + Ger = — a1)(y’ — ae) 
Y + + = Vly’ — Bi)(y’ — Be). 


Their resultant R is 
R = (a1 — B1)(a1 — B2)(@2 — B1)(a2 — Be) = — 


where 


H= + WorPyy 


The condition that the tangents named be not parallel is R + 0, or since 
Ag = 0,H +0. Asin IIIa, 2 and 3, the nature of the extreme is given by 
the sign of gy, if gy, + 0; by the sign of gz if gy, = 0; not both ¢,, and 
¢y vanish, and they have the same sign if neither vanishes. If H = 0 
the two tangents are parallel and the discussion fails. 

3. AY > 0. The curve C has two real branches intersecting at P with 
distinct tangents. It is evident geometrically that wo is always an extreme 
of ¢ for one branch of C, and an extreme of the same kind for both branches 
if the tangent to neither branch is parallel to the asymptotic tangent to S 
at @. Considered analytically, we have R = H?, since Ag = 0, and a; = az. 
If H + 0, we have for the two branches of C, since y,’’ is finite, d’u/dz* 
equal to gy(81 — ai)? and ¢y,(B2 — a1)”, gy +0; the two extremes are 
alike, maximum or minimum as ¢,, is negative or positive. If H = 0 the 
tangent to one branch of C is parallel to the asymptotic tangent at Q, and 

yr this branch the discussion fails; for the other branch the nature of the 
extreme is given as before by the sign of gy). 

IIIc. Ag > 0. The total curvature of S is negative at Q; there are two 
distinct asymptotic tangents. 


i 
H 

4 
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1. AY <0. There is no extreme of ¢ at wp. 

2. AY = 0. It is evident geometrically that sufficient conditions that 
Uo be a relative extreme of ¢ are that P is not an isolated point of C and that 
the single tangent to C at P is not parallel to either asymptotic tangent at 
Q, that is H +0. If H = 0 the discussion fails. 

We determine the nature of the extreme when H + 0 as follows: In 
the neighborhood of 2, yo, we have 


g= Ut (Ay)*Gyy + 2ArAygzy + 


for values of the partial derivatives at a point near 20, yo. Since the deriva- 
tives are continuous by hypothesis wo is a maximum or minimum as 


Puy + Gey + Pex = 01)(B1— 02), = Bi = Br, 
is negative or positive. Now it is easily proved that 
H = (a1 — Bi)(a2 — Be) + (a1 — — Bi)}, 
and, since in the case before us 8; = Bo, we have 
H = Gy — Bi)(a2 — 1). 


Then wp is a maximum or minimum as y,,H is negative or positive. The 
value of d?u/dx? cannot be used in this discussion since in general, when 
Ay = 0, y,’”’ becomes infinite. 
3. AY > 0. For both branches of C y,’’ is finite if Y,, + 0, and, if 
Yy + 0, 
du 


dz? = Gy(y’ — a1)(y’ — ae). 


If R < 0, the values of this second derivative for the two branches of C, 
— a1)(B1 — a2), Pryy(B2 — a1)(B2 — ae), 


have opposite signs, and wp is an extreme of ¢ of opposite kinds for the two 
branches of C. If R > 0, the two values have the same sign and wp is an 
extreme of ¢ of the same kind for the two branches of C. If R = 0, one or 
both tangents at P are parallel to respectively one or both asymptotic 
tangents at Q. Along such a branch of C the discussion fails. If both 
tangents at P are parallel to asymptotic tangents at Q we must have 


Grr Pry _ Py 


Yor Vey Vw 
We consider the nature of the extremes of g at uw. Suppose first R > 
there are two extremes of the same kind. If yg, and W,, are both different 
from zero these extremes are maxima or minima as the following expression, 


¢ 
fi 
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not zero, is negative or positive: 
9 
{ (a1 — B1)(a2 — Bi) + (a1 — B2)(a2 — B2)} = U4, 
yy 
= + 2¢yAy. 


The nature of the extremes is then given by the sign of Uy. It will 
appear in § 4 that the restriction ¢,,~,, + 0 is not essential. If R = 0, 
a; = Bi, a + Bo, the criterion for maximum or minimum on the branch 
of C for which y’ = z is again the sign of U4, unless gy, = Py, = Us, = 0, 
when the distinction will be given by the sign of 


U,’ Worl + 0. 


We shall show in § 4 that U, and U,’ have the same sign when neither 
vanishes. 

In three other cases, previously considered, the nature of the extreme 
or extremes is given by the sign of U, or of U,’ when the former vanishes. 
These all come under III, 9, = g, = vz = Yy = 0; they are b, Ag = 0, 
2, Ay = 0; 6b, 3, AY >0; c, Ag > 0, 2, AY =0. In the discussion of 
b, 2, H + 0, it was shown that wp is a relative maximum or minimum as ¢y,, 
supposed different from zero, is negative or positive; in this case Uy, = Wy,H 
and this expression has the same sign as ¢,, if ¥,, does not vanish, for 


= — Bi)? > 0. 


In 6, 3 the criterion was again the sign of ¢,,, which has the same sign as 
U,, when these are both different from zero, whether or not H vanishes; 
if either H or y,, vanishes U, = g,Ay; if Y,, + 0, we have from the preced- 
ing paragraphs, 

{ (a1 — Bi)? + (a1 — Be)?} = 


In c, 2 it was proved that if H + 0 the extreme is maximum or minimum 
as W,,H is positive or negative; since Ay = 0, this is exactly U,. If H = 0, 
U, = U,’ = 0, and the discussion fails. If Ag = AY = R= 0 we have 
the excluded case, 


Par _ Puy _ Puy, 


If Ag > 0, AY = R = 0 the discussion fails. In all other cases where 
Ag 2 0, Ay 2 0, gz 0, 


gy has at uw a maximum or a minimum as U, is negative or positive; if 
U, = 0 the criterion is the sign of U,’ + 0. In certain cases ¢ has at wu 
two like extremes whose nature is given by the sign of U, or U,’, namely 
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Ag = 0, AY > 0, R> 0. In other cases ¢ has at wp one extreme whose 
nature is given by the sign of U, of U,’, and may have a second extreme of 
either kind or no other extreme, namely Ag = 0, Ay > 0, R = 0. 

§ 3. Reciprocity in the Exceptional Cases.—We now consider problems 
A and B together under the same assumed conditions with the purpose of 
determining, in so far as it is possible with the use of the first and second 
derivatives of ¢ and y, when extremes in the two problems are like or unlike. 

It is clear that under assumed conditions for example in case I of § 2 the 
solution of B is given by interchanging ¢ and y in II. 

In I, yz = Wy = 0, ¢- and gy, are not both zero. Ia, AY <0. In A 
there is no extreme. For B the solution is given by Ila, exchanging ¢ 
and y; there is always an extreme. 

Ib. AY = 0. In A if P is not an isolated point of C a sufficient condition 
for an extreme is Dyg + 0. For B the solution is that of IIb; a sufficient 
condition is Dyg + 0. The discussion fails for comparison of the nature 
of the extremes in A and B. The sufficient condition common to the two 
problems is, geometrically expressed, that the two curves, ¢ = uo and 
y = %, are not tangent at P. 

Ic. AY > 0. In A it is necessary for an extreme that Dyg = 0. For 
an extreme in B a sufficient condition, by IIe, is Dyg + 0. 

The solutions of A and B with the assumed conditions of case II are 
given by exchanging ¢ and y in the preceding paragraphs. 

In I, g = gy = te = Wy = 0. 

IIIa. Ag < 0. 1. AW < 0. There is no extreme in A or B. 

2. AY = 0. There is always an extreme in A unless P is an isolated 
point of C. The solution of B is given by IIIb, 1; there is no extreme. 

3. AY >0. There is an extreme for any path in A. In B, by IIIe, 1, 
there is no extreme. 

IIIb. Ag = 0. 2. AY = 0. In both A and B, if P is not an isolated 
point of either curve, ¢ = wp or C, a sufficient condition for an extreme is 
H +0. We prove that the extremes are like or unlike as H is positive or 
negative: Since Ag = Ay = 0, ai = ae and 6; = B2; since R = H? + 0, 
a, + B, and H = gyy(a1 — Bi). Then ¢g,, and y,, have like or unlike 
signs as H is positive or negative. 

3. Ay > 0. In A there is always one extreme, and if H + 0 two like 
extremes. In B, if P is not an isolated point of the curve ¢ = wo, a suffi- 
cient condition, by IIIc, 2,is H + 0. If H + 0 the extremes of A and B 
are like or unlike as ¢,, and H¢y,, have like or unlike signs, that is as H is 
positive or negative. 

IIIb. Ag > 0. 3.4y¥ > 0. If R < 0 there are, as has been shown in the 
preceding section, two unlike extremes in problem 4, and similarly two 
unlike extremes in B. 
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For R = 0, we have, with the notation of § 2, a; + a2 since Ag > 0, 
and 6; + Be since Ay > 0, and a; = f). If also az = fy then 


Por _ Pry _ 
Ven Vy Vy. 

and further discussion fails in both A and B. If these equations do not 
hold a2 + 2 and there is certainly one extreme for both problems. For 
A and B respectively, if 9,,,, + 0, 


2 
= — a1)(B2 — az), = Wy(a2 — B1)(a2 — B32), 


and since, for this case, 
H = — B2)(a2 — Bi), 


the extremes known to exist in A and B are like or unlike as H is positive or 
negative. 

For R > 0 there are two like extremes in A, similarly two like extremes 
in B. We shall show again that the.extremes of A are like or unlike those 
of B as H is positive or negative. Supposing as before ¢,,/,, + 0, the two 
values of d?u/dx? for the two branches of C in A have the same sign and 
are, as given in § 2, 


Pyy(B1 — a1)(B1 — a), Pyy(B2 — o1)(B2 — a2). 
The two values of d?v/dz? in B are similarly 
— Bi)(a1 — Be), — B1)(a2 — Be) 


and have the same sign. The signs of the two products, 


— B:)(a2 — Be), PyyWyy(@1 — Be)(a@2 — Bi), 


are alike and the same as the sign of H since the latter is one half their sum. 
Then the extremes of A and B are like or unlike as H is positive or negative. 

We note that in all of the exceptional cases of § 2, where it is possible to 
determine the nature of the extremes, supposed alike when there are two, 
in both A and B by the use of the first and second derivatives of ¢ and y, 
the extremes in the two problems are the same or different as H is positive 
are negative. These all come under III, ¢, = g, = ¥z = ¥, = 0, and 
include all cases where neither discriminant, Ag or Ay, nor the resultant R 
is negative, not all three vanish, and H + 0. 

§ 4. Invariant Properties.—The various analytic conditions derived in 
the preceding sections must be independent of the geometrical representa- 
tion of the problem and are accordingly invariant for any change of variables, 


be 
& 
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not singular at the point considered. It is of interest to consider directly 
the invariant properties of these conditions. It is also necessary to demon- 
strate the unessential character of assumptions as to non-vanishing of 
certain derivatives made in some of the discussions. 

Suppose 2 and y are replaced in ¢ and y by two real functions of two 
new real variables, x’ and y’._ Concerning these functions we assume that 
both have for the point considered finite first and second partial deriva- 
tives; further, that at the same point the transformation is not singular, 
that is 

6 = — + 0. 


For the new variables the various conditions will be obtained by replacing 
the derivatives of ¢ and y with respect to x and y in the conditions estab- 
lished by the corresponding derivatives of the same functions with respect 
toz’ andy’. We have 


Px! = Pate’ + Py! = + 

= Prats! + + + + 

Px'y’ = F Pry(Le'Yy’ + Ly’Yr') + + + 
Py'y’ = + + + + 


with similar equations for y. 
Consider the quadratic form, 


Ap? + 2Bpq + Cq? = A’p” + + C’q”, 
p=ap'+ by’, q=cp'+dy, 
A’ = Aa’ + 2Bac+ Ce’, 
B’ = Aab + B(ad + be) + Ced. 
C’ = Ab? + 2Bbd + 
If we write 
a= 2,’, b= C= Yz', d = yy’, 


we may consider ¢, and ¢, as variables contragredient to p and q.* In the 
following discussion of the invariant properties of the conditions of the 
preceding sections it is to be remembered that the values of all derivatives 
of y, y, x, and y considered are the values at P (2x0, yo). In the first place 
the classification of § 1 and in § 2, cases I, II, III is invariant, for, since 
5 + 0, g2’ = gy,’ = 0 if and only if ¢, = g, = 0. Next, if 9, and g, are 
not both zero, and if yz and y, are not both zero it is clear that x’ and y’ 
may be chosen so that ¢,’ and y,’ are both different from zero. The 
necessary condition of § 1, J-= 0, is invariant for J is a contravariant* of 
weight one. The condition of the same section, sufficient for an extreme 


* Bécher, “Introduction to Higher Algebra,” p. 109. 
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and determining the nature of the extreme, we consider later. When the 
sufficient condition is satisfied for both A and B the condition that the 
extremes are like or unlike is that g, and y, have unlike or like signs. If 
y, and y,’ are both different from zero this condition is invariant if J = 0, 
for 

bez + doy _ vy 


bztdy 
We may say that ¢,/¥, is an absolute conditional invariant, with condition 
J = 0. 

In § 2 the further subdivision of the problem depends on the signs of 
Ag and Ay, Ag = ¢n,/? — GrrGy,. It is evident from the values given 
above for the second derivatives of ¢ with respect to x’ and y’ that Ag is 
not in general invariant, but with the conditions, g, = ¢, = 0, the trans- 
formation is exactly the transformation of the quadratic form in p and q, 
and since the discriminant is an invariant of weight two in the algebraic 
theory,* we have conditionally A’g = &Ag. The sign of Ag and the 
vanishing of Ag, and similarly Ay, are therefore conditional invariants: 
If = = 0, the expression, 


= Vyy = + Wer, 


is invariant in sign, for it is in the algebraic theory the form adjoint to 
Warp? + Weypq + Wyyg’, and therefore a contravariant of weight two.t 

We now consider the invariance of sign of g,(y,’”’ — y4’’), supposed 
not zero, in §1. We may evidently write 


Wy y 
where we suppose ¢, and y, both different from zero, and J = 0. Under 
these conditions the sign of the expression is evidently unchanged by trans- 
formation if g,’ and y,’ are also both different from zero, since ¢,/y, is an 
absolute invariant, and since the terms of D,’y and D,’y, which contain 
the second derivatives of x and y with respect to 2’ and y’, cancel. 

In case III we have the conditions, ¢, = yg, = ¥z = ¥, = 0, so that the 
transformation is equivalent algebraically to that of two quadratic forms. 
We have supposed in that discussion that the following equations do not 
hold: 


Wer Vay Yuy 


* Bécher, 1. c., p. 129. 
t Bécher, 1. c., p. 159. 
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which include the special cases also excluded, 


= Pry = Cyy = 0, or Ves = Wry = Vuy = 0. 


That this condition is invariant when the four first derivatives vanish is 
obvious from the values given above for the second derivatives of ¢ with 
respect to x’ and y’. From the same expressions it follows that under the 
assumed conditions it is possible so to choose 2’ and y’ that neither ¢,, 
nor Y,, vanishes. In the same case we have stated that the nature of the 
extreme is in certain subcases dependent on the sign of gy, or Wyy. We have 


Py'y’ = + + 


and the signs of ¢,, and ¢,’,’ are the same when this form is definite, 
Ag < 0, or, if neither vanishes, when the form is singular, Ag = 0. Still 
considering the same case we have from the algebraic theory* the facts that 
H and R are invariants of weight two and four respectively, and are therefore 
unchanged in sign by transformation. | 
It remains only to consider the invariance of the sign of 


Us Vy + 2¢,Ay, 


the criterion for maximum or minimum in III, c, 3, Ag > 0, Ay > 0; 
R20. We recall the fact that both H and Ay are invariants of weight 
two, and note that H does not vanish since R = H? — 4AgAyY = 0. Nowif 
¢y and yy, are both different from zero it has been proved that the sign of 
U, is the same as that of 


{ (a1 — Bi)(a2 — Bi) + (a1 — Be)(a2 — B2)}, 
and U, + 0. Evidently if one but not both of the derivatives, gy, and Py, 
vanishes U, + 0. If also ¢,: and yz; are both different from zero the sign of 


U4! = Wer + 


is the same as that of 


Hava 
Bx By a, a2 Be 
Puy (a1 — Bi)(a2 — Bi) — — Be) 
B? 


and consequently in the cases before us, R 2 0, the same as that of U,. 
Any given non-singular transformation, 


p=ap'+ by, q=cp'+dy, 
* Bocher, |. c., pp. 166, 236. 
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may be regarded as the transformation, 6 = 1, of the transformations @, 
p= a(0)p’+ b(0)q’, g=c(O)p’+ d(A)q’, 


where the four functions of 6, a, b, c, d, are single-valued functions, continuous 
for 0 S 6 = 1, and subject to the conditions: 


a(0)=1, a(1)=a; 6(0)=0, =); c(0) = 0, c(1) =¢; 
d(0)= 1, d(1) =d. 


Evidently @ = 0 gives the identical transformation. If for the given 
transformation ad — be is positive the four functions of 6 may be chosen 
so that in the interval named 


a(6)d(6) — b(8)c(@) 


does not vanish, and consequently no transformation @ is singular. For the 
same interval U, is a continuous function of 6. If R > 0 not both ¢,, 
and y,, vanish for any 0, and U, has the same sign for 6 = 0 and @ = 1. 
If for the given transformation ad — bc is negative let us suppose, as we may 
without loss of generality, that originally no one of the four derivatives 
Wer» is zero. Let the transformation with 6 = — 1, 


be first applied; U, is transformed without change of sign to U,’, and may 
then be proved, as in the case of positive determinant, to be unchanged in 
sign by the given transformation. 

The case R= 0, a: = 81, a2 + Bo, needs further consideration, for 
U, = 0 if gy = vy = 0. If however Us, = 0 necessarily U,’ + 0, for 
the possibility, 


Prz = Pyy = Ver = 0, 


is excluded. It remains to show that U, and U,’ have the same sign when 
neither is zero. Suppose, as before, that for 6 = 0 none of the four deriva- 
tives vanishes; then U, and U,’ have the same sign. When @ increases 
from 0 to 1 U, and U,’ have continually the same sign unless one, say U, | 
vanishes for 6 = 6;, when we must have ¢,, = Y,, = 0. It is now con- 
ceivable that, for further increase of 6, U, should change sign if either ¢,, 
or Yy, remains zero while the other is not zero. We prove that, with our 
assumptions, this cannot occur. Suppose, for ? = 01, gyy = Wy, = 0; that, 
for all values of 6 such that @ — 6; is positive and less than some positive 
number ¢, gy, = 0, yy +0. From R = 0, we have 


Vyy 4 + 40 War = 0, 
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for any value of 6 between 6; and 6; + «. This equation is satisfied by 1. 
rr = Wzr = 0. This is inadmissible, since we cannot have, for 41, G22 = Yyy 
= Wor = Wy = 0; 2. Ger = Ory = 0. This case is also excluded, for we 
cannot have, at 01, G22 = Gry = Sy = 0; 3. if ger + 0, 
ry 

Since all the derivatives are continuous functions of @ and since, for 41, 
Wyy = 0 we should have in this case for 6,, and consequently for all trans- 
formations, again the excluded case 


Consequently our hypothesis is inadmissible, and U, is always of the same 
sign when it is different from zero. In the case III, b, 2, 9, = gy, = yz 
= y, = 0, Ag = Ay = 0; H +0, wehave Uy =yy,H. It has been proved 
that y,, cannot change sign though it may vanish, but not both y,, and zz 
vanish. In b, 3, Ag = 0, AY > 0, we have, if H = 0, Us, = gyAy and, 
since gy, may vanish but cannot change sign, the same is true of U,; not 
both gy and gz, vanish. If H + 0 the proof given above for c, 3, R > 0, 
applies, and U, can neither vanish nor change sign. For c, 2, Ag > 0» 
Ay = 0; H + 0, we have, as in b, 2, Us = YH. If H = 0 the discussion 
fails. 

We remark that our discussion proves that U, and U, have the same or 
different signs when neither is zero in the cases considered as H is positive 
or negative. 

It is of interest to note that our discussion proves the existence of certain 
absolute conditional invariants or conditional differential parameters, for 
example 


the first with the conditions, ¢, = ¢, = 0, the second and third with the 
conditions 9g; = = = = 0. 

§ 5. Examples.—We give in this section simple examples illustrating 
each case of the theory developed in the preceding sections, also examples 
showing, in each case where the discussion is said to fail, that there may 
exist a maximum or minimum or no extreme. We shall illustrate the 
reciprocity in problems A and B, and the invariance of the conditions in 
some of the examples by interchanging the variables x and y. All of the 
examples are so chosen that x = yo = Uo = % = 0. 
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For the general case of § 1, ¢; and ¢, not both zero, ¥, and y, not both 
zero. 

1. 2. 
We have gy, = ¥y = 1; J=0; ys” = 0, y,” = 2. The curves, y = 0 and 
y = 2’, are tangent at (0,0) but do not osculate. In A we have a relative 
minimum, in B a relative maximum; ¢, and y, have the same sign. 


2. y, y=y— 2, #, y+ 2. 


For each y in this example the curves, y = 0 and x = 0, osculate at the 
origin; for the three choices ¢ has respectively no extreme, a minimum, a 
maximum. 

For § 2, case I, ¢, and ¢, not both zero, yz = py = 0. 

a. Ay < 0. 

3. g= y, 

Evidently no real branch of C, y = 0, passes through the origin. 

b. AY = 0. The point P (0, 0) is a cusp, an osculating point, or an 
isolated point of C. 

For the first two y, the origin is a cusp of C, and Dyg = 2. For the first 


y there is a minimum, for the second a maximum of ¢. For the third choice 
of y the origin is an isolated point of C and ¢ has no relative extreme. 


5. g=y, y = 7? — y’, (y a”)? x5, (y+ a’)? 


For each choice of y the origin is a cusp of C, and Dye = 0. For the three 
choices ¢ has respectively no extreme, a minimum, a maximum. 


6. g=y ez. 
The origin is an osculating point of C, and Do = 0. For a path on C 
with continuously turning tangent ¢ has at the origin no extreme, a mini- 


mum, or a maximum, depending on the path. 
ce. AY > 0. 


For each value of y the origin is a double point of C with two distinct tan- 
gents, and Dyg = 0. If it be desired that C be irreducible y? may be added 
to each value of y without affecting the results. For the three choices of y 
the value zero is respectively no extreme, a minimum, a maximum of ¢ 
for that branch of C whose tangent is not parallel to the horizontal tangent 


to S at Q. 
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For § 2, case II, g, = g, = 0, ¥z and yy, not both zero. 

a. Ag < 0. 

8. y = y. 
The two choices of ¢ give D,y the values 2 and — 2 respectively. For the 
first ¢ has a minimum, for the second a maximum. 

b. Ag = 0. 

9. g=xr— 7, yf, y=y. 
Both choices of ¢ give Dyy = 2, and in both cases g has a minimum. It is 
evident on comparing this example with example 4 that we have no basis 
for comparing the extremes of problems A and B with the hypotheses of Ib. 
Similarly if we exchange ¢ and y as given in examples 5 and 6 we have for 
each D,y = 0, from 5 no extreme and two minima, from 6 a maximum. 

c. Ag > 0. 

10. yz’, y= 2, y. 


The two values of ¥ give D,y equal to — 2 and 2 respectively and give 9 
maximum and minimum values respectively. For problem B, correspond- 
ing to Ic there is no extreme. If we interchange ¢ and y as given in example 
7 we have D, and Y= 0 and have a maximum in the first case, no extreme in 
the second and third. There is evidently no possibility of comparing the 
nature of the extremes in A and B in this case furnished by our discussion. 

For § 2, case III, 9, = vy = vz = ty = 0. 

a. Ag <0. 1. Ay <0. 


11. g=xr?+y, y = 227+ y’. 
There is no relative extreme in either A or B. 

2. Ay = 0. 

12. att 


Let y be given any of the values assigned in examples 4, 5, 6. For all 
choices ¢ has a minimum and gy, = grr = 2. 
3. Ay > 0. 
12. 
Let y have the values of 7. There is a minimum of ¢ for each branch of C. 
b. Ap = 0. 1. AY <0. 


14. g=¥, 
Clearly ¢ has no relative extreme at the origin. 

2. Ay = 0. 

15. g= y= 2’. 


We find yy = 2, H = 4, U, = 0, U,’ = 8, U, = 8, U,’ = 0. 


= 
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Problems A and B have like extremes, both minima. 


16. g=-—y’, y = 2’. 


Cy = — 2, H= —4, U, = 0, U,'’ = — 8, U, = 8, U,’ = 0. 


Problems A and B have unlike extremes, a maximum in A, a minimum 
in B. 

The case, Ag = AY = R= H = 0, is excluded, since, as previously 
stated, we have, when these all vanish, 


Vor Vy Vy 


The following example shows that in this case g may have no extreme, a 
minimum or a maximum : 


17 g=y¥ y=y— 2’, 
3. Ay > 0 
18. g=y, yp=or-y¥ 
We find 
Cyy = 2, Vy = — 2, Ay = 4, H =4, 
U, 8, U,’ 8, U, 8, U,’ 0 
Problems A and B have like extremes, both minima. 
19. y’, y= y 
We have 
Cy = 2, Vy = 2, Ay= 4 H =-4, 
U, = 8, U,’ = 8, U, = — 8, U,’ = 0. 


Problems A and B have unlike extremes, minima in A, a maximum 
in B. 

20. y— y = zy. 
We have for all choices of ¢ 


Cy = 2, Yy = 0, Ay = 1, H =0, 
U, = 4, U,’ 0, Uy, 0, U,’ 0. 


I 


Problem A has for the three choices of ¢ respectively a minimum and 
no extreme, two minima, a maximum and a minimum. In each case the 
nature of the known extreme, a minimum, is given by the positive sign of 
U,. The discussion fails for B. : 

c. Ag >0. 1. AY <0. There is no extreme in problem A. 
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2. Ay = 0. 

21. g=P-y, y=y¥y, 
We find H = 4, U, = 8, U,’ = 0. There are a minimum in problem 4, 
two minima in B. See example 18. 


22. g=a+2, xy — y=y, y— 2’. 
We have H = U, = U,' = 0. For the three choices of ¢ problem A has, 
for a path C with continuously turning tangent, respectively no extreme, 
a minimum, anda maximum. The discussion fails to give any information. 
See the reference to problem B in example 20. 

3. Ay > 0. 
23. = 2a? — + 2y’. 


We find Ag = 4, Ay = 9, H = 0, R = — 144. In problem A the values 
of ¢ for the two branches of C are 3y? and — 32”, a minimum and a maximum 
respectively. In B the two values of y are 9y? and — 2”, a minimum and a 
maximum. 


24, 

We find 
Ag = 4, Ay = 16, H = — 20, R= 14 
U, 96, U,’ 24, U, 24, U,’ = — 24, 


In problem A there are two like'extremes, since R > 0, minima since U, > 0, 
in B two like extremes, maxima since U, < 0, unlike those of A since 
H<0. 
25. y = — 2’. 
We have 

Ag = 4, Ay = 16, H = 20, R 144, 

U, = 96, U,' = 24, U, = 24, U,’ = 24. 


There are two like extremes, both minima, in both A and B, of the same 
kind since H > 0. 


We have 
Ag = 4, Ay = 1, H = —4, R= 0, 
U, = U,’ = 4, U, = U,’ = — 8. 


For that branch of C whose tangent is not parallel to an asymptotic tangent 
to S at Q, that is for y = 0, ¢ has a minimum value, since U, > 0; for the 
other branch, x — y = 0, the three choices of ¢ give respectively no extreme, 
a minimum, and a maximum. In problem B the branch of the curve, 
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gy = 0, whose tangent, x + y = 0, is not parallel to an asymptotic tangent 
to the surface, y = 0, at the origin, gives y a maximum value, since U, < 0, 
unlike the extreme in A since H < 0. 


27. x’, 
We have 
Ag = 4, Ay = 1, 
U, = U,' = 4, 


In this case, since H > 0, the extremes, different in A and B in example 26, 
are the same, minima in both problems. a 
March 28, 1920. 
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y= y — xy. 

| Red, 
U, = U,,’ = 8. 
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